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COMBUSTION OF SOLID CARBON RODS IN ZERO AND NORMAL GRAVITY 


Charles M. Spuckler 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

I. INTRODUCTION 

In recent years, a new emphasis has been placed on research in- 
volving the production of energy from new or existing sources. One of 
in areas of interest is solid combustion oriented toward an under- 
standing of the true mechanisms involved in the combustion of coal. 

The chemical and physical processes involved in the combustion of a 
solid depend in part on the nature of the solid. Many solids burn in 
a manner similar to liquid fuels, in that the surface of the solid 
sublimates and the burning takes place away from the surface in a 
flame zone surrounding the solid. Other solids burn with a surface 
reaction. If the surface reaction results in complete combustion, no 
flame zone will appear but the surface will glow. If combustable 
gases are formed by the surface react ion, they can burn in a gas phase 
reaction producing a flame zone away from the surface. 
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Combustion of carbonaceous solids, such as coal, is complex be- 
cause it takes place in two stages. The coal first pyrolyzes on being 
heated, with the volatile gases being driven from the solid through 
pores and cracks to the outside environment where they mix with oxygen 
and burn via a gas phase reaction. Auer all volatiles have been 
driven off, mainly solid carbon remains, which then burns via re- 
actions on the surface and in the pores. Ouring this latter process, 
various chemical reactions can occur. The carbon surface reacts with 
oxygen to form either carbon monoxide, carbon dioxide, or a combin- 
ation of both. If carbon monoxide is formed at the surface, a gas 
phase reaction of the carbon monoxide with oxygen to form carbon 
dioxide is possible. 

To gain a better understanding of the mechanism of coal combus- 
tion, it was decided that the process by which carbon ignites and 
burns should be examined. If the carbon combustion process can be 
understood, the results should aid in the analysis of the larger prob- 
lem of coal combustion. 

In this study, zero and normal gravity experiments were conducted 
to determine the correct combustion model for carbon burning supported 
only by its own combustion. Species concentration profiles were meas- 
ured to determine the mechanism controlling the process and the pri- 
mary products of combustion. The zero and normal gravity data were 
then compared to ascertain the effects of free convection. 


II. LITERATURE SURVEY 


A. General Discussion 

The term heterogeneous combustion is applied to any combustion 
process in which a fuel in the condensed phase is burned. This defi- 
nition covers the following combustion processes (1) tne fuel is vap- 
orized at the surface and the chemical reaction occurs in the gaseous 
phase; and (2) the chemical reaction occurs on the surface of the con- 
densed phase with possible secondary reactions occurring in the gas- 
eous phase. Liquids and some solid fuels fall into the first cate- 
gory. Tho fuel is vaporized at the surface and moves outward, meeting 
oxygen which is moving inward toward the surface, with the combustion 
occurring in the gas phase away from the surface. For ease of solving 
this problem analytically, the chemical reactions in the gas phase are 
generally assumed to occur infinitely fast, so that analytically the 
combustion occurs in a very thin flame sheet. The products of com- 
bustion move both inward toward the fuel surface and outward from the 
flame zone. Some of the heat produced by the reaction is transferred 
to the surface where it vaporizes more fuel, while the rest of the 
heat is transferred away. A schematic of the process is shown in Fig- 
ure 1. 

For solids which fall into the second category, the combustion 
process consists of oxygen first being transferred to the surface. 
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where it is chemisorbed and where it reacts with the fuel. The pro- 
ducts are then desorbed and move away from the surface. If the chemi- 
sorption, the surface reaction, and the desorption processes are lump- 
ed into one overall chemical process, solid combustion can be said to 
consist of a mass transfer process and a chemical process. 

Carbon is one of the solids that burns with a surface reaction 
and, therefore, falls into the second category. Various models for 
carbon combustion have been proposed. The products of combustion and 
the mechanism of the combustion process depend upon the carbon surface 
temperature. 

B. Theoretical Analysis 

The various models proposed for carbon combustion can be classi- 
fied as a low temperature model, a mid-temperature model, and a high 
temperature model. Most of the models found in the literature are for 
a quiesent atmosphere and neglect buoyancy effects. 

At low surface temperatures, below approximately 500* C, carbon 
is oxidized at the surface to carbon dioxide through the reaction 
C + 0^ » CO^. This model along with the assumed concentration and 
temperature profiles is presented in Figure 2. Fendell [1J presents a 
general solution for this model along with limiting solutions for mass 
transfer control and chemical process control. Also, solutions for 
equilibrium and nonequilibrium reactions at the surface are presented. 

At moderate temperatures, between 500° C and 1100* C, a complex 
reaction takes place in which both carbon dioxide and carbon monoxide 
are formed at the surface through two surface reactions, C + Op » 


OXYGEN 


U 

CARBON OIOXIDC 


SURFACL REACTION: C ♦ - CO 2 



Figure 2. - Low temperature model. (T s < 58^ C. ) 



7 


C0 2 and 2C ♦ 0 2 * 2C0. There is no oxidation of carbon monoxide 
in the gas phase. Ther Is no known theoretical analysis of this 
model, which is illustrated in Figure 3. 

For high surface temperatures, above 1100* C, four high temper- 
ature models have been proposed. Ubhayahar and Williams [2, 3] pro- 
posed a model in which oxygen reacts with carbon at the surface to 
form carbon monoxide through the reaction 2C ♦ 0^ » 2C0, with no 
further reaction of carbon monoxide occurring in the gas phase. This 
model is based on extinction experiments in which carbon particles 
were burned in oxygen, and oxygen and nitrogen mixtures. In their 
theoretical analysis, atomic oxygen in equilibrium with molecular oxy- 
gen or just molecular oxygen reacting at the surface were considered. 
This model, neglecting the dissociation of oxygen, is presented in 
Figure 4a. 

Fendell [1] proposed a complex high temperature model in which two 
surface reactions take place. Carbon is oxidized to carbon monoxide 
through the reaction 2 C + 0^ » 2C0, and ambient carbon dioxide is 
reduced to carbon monoxide through the reaction CO^ + c * 2C0. 

There is no gas phase oxidation of carbon monoxide. Fendell presents 
the general solution along with the limiting solutions of mass trans- 
fer control and chemical process control. This model is presented in 
Figure 4b. A third high temperature model has been proposed. In this 
model carbon monoxide is formed at the surface through two surface 
reactions, C + » 2C0 and 2C + 0 2 » 2C0. The carbon monoxide 

is then burned to carbon dioxide through the homogeneous gas phase 
reaction, 2C0 + 0,, » 2C0 2 . This model is presented in Figure 4c. 

Hugo, Wicke, and Wurzbacher [4] obtained an approximate solution for a 







CARBON 


GAS PHASE REACTION 



SURFACE REACTION; C ♦ CO? * 2CO? 
GAS PH&SE REACTION: 2CO ♦ 0? * 2CO? 
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GAS PHASE 
REACTION RATE 


PROFILES FOR 
FINITE RATE 
GAS PHASE 
REACTION 


0 r - 

DISTANCE FROM SURFACE 

Figure 4. - Concluded. 
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cylindrical carbon rod burning according to this model. The gas phase 
burning was assumed to be a first order reaction with respect to car- 
bon monoxide. The two simplifying assumptions that made the solution 
approximate were that the carbon monoxide concentration in the bound- 
ary layer was small and that the boundary layer was thin so that an 
average boundary layer thickness could be used. The theoretial carbon 
monoxide concentration profiles showed good agreement with the experi- 
mental data obtained by Wicke and Wurzbacher [5]* who burned vertical 
carbon rods in oxygen. Caram and Amundson [6] obtained numerical so- 
lutions for carbon slabs and spheres burning according to this model. 
They assumed that the gas phase reaction was first order with respect 
to carbon monoxide, one-half order with respect to oxygen and one-half 
order with respect to water vapor. The ambient concentration of water 
vapor, which was considered small enough so that other possible re- 
actions could be neglected, was considered a parameter of the model. 
They also determined feasible surface temperatures. The feasible sur- 
face temperatures were those that allowed the mass fractions to be 
equal to or greater than zero and less than or equal to one. 

A fourth high temperature model has been proposed. In this model, 
carbon dioxide is reduced to carbon monoxide at the surface through 
the reaction CC >2 + C > 2C0. The carbon monoxide moves outward, 
meeting oxygen, and is burned to carbon dioxide through the gas phase 

reaction 2C0 + » 2CC>2. Analytically, the gas phase reaction 

can be considered to occur either with an infinitely fast reaction 
rate so that a very thin flame sheet forms, or with a finite rate so 
that a flame zone appears. This model is presented in Figure 4d. 
Spalding [7] presents a solution for this model using the assumption 
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that thermodynamic equilibirum exists at the surface and throughout 
the gas phase. Fendell [1] also considered this model and presents 
solutions for two different cases: (1) mass transfer controlled, 

equilibrium, irreversible surface and gas phase reactions; and (2) a 
finite rate surface reaction and a mass transfer controlled gas phase 
reaction. Van der Held [8] presents a l»ion for this model for a 
one dimensional caroon surface burning flowing gas. He consider- 
ed the case of an equilibrium surface carbon dioxide reduction re- 
action with a homogeneous carbon monoxide oxidation reaction of either 
2/3 order, second order, or third order. 

The burning rate equations that are generally used for carbon com- 
bustion are derived by first solving the continuity equation for the 
flux of oxygen at the carbon surface. The partial pressure of oxygen 
at the surface is eliminated from the oxygen flux equation using an 
Arrhenius rate equation. This results in a burning rate equation. 
Burning rate equations that have been obtained for carbon spheres and 
cylinders are presented in Table 1. The definition of the variables 
in the equations can be found in the Nomenclature Section. 

C. Controlling Mechanism 

The general burning equations for carbon, such as those in Table 
1, are of the form: Combustion Rate * 1 /(Rq + R c ) where Rp refers 
to the mass transfer resistance and R c refers to the chemical pro- 
cess resistance. Therefore, the combustion process can be considered 
to consist of two processes, a chemical process and a mass transfer 
process. The mass transfer process consists of convective mass trans- 
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Table 1* - C .nued 
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Table 1* - Continued 
spheres 
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fer and molecular diffusion, and the chemical process consists of 
chemisorption, the surface reaction, and desorption. The overall com- 
bustion rate of a solid particle will thus be determined by the slower 
of the two processes, that is, either the chemical process or the mass 
transfer process will control the burning rate. 

Various parameters affect both the chemical process rate and the 
mass transfer rate. The chemical process rate is a function of the 
material, the surface temperature, and, if the reaction is not zero 
order, the oxygen partial pressure. The mass transfer rate is a func- 
tion of ambient gas velocity and temperature, diameter and temperature 
of the particle, and oxygen partial pressure. As the particle di- 
ameter decreases, the mass transfer resistance decreases and a criti- 
cal diameter is reached where the chemical process controls the com- 
bustion rate. Similarly, for a given diameter particle, as the gas 
velocity increases the mass transfer resistance decreases, and at some 
critical velocity the chemical process will control the combustion 
rate. Mulcahy and Smith [16] state that for a fixed size particle, 
increasing the temperature increases the chemical reaction rate and 
therefore the chemical process rate, and at some critical temperature 
the ma .s transfer will control the rate of combustion. Ubhayakar and 
Williams [2, 3] conducted extinction experiments and stated that ex- 
tinction was caused by a change from mass transfer control to chemical 
process control. They found that as the partial pressure of oxygen 
was increased, the particles burned to smaller diameters. Therefore, 
it can be concluded that as the partial pressure of oxygen increases, 
the critical diameter at which the chemical process controls the com- 
bustion rate decreases. 


18 


The temperature and particle size, at which a change from chemical 
process control to mass transfer control has been found to occur, are 
presented in Table 2. From this table, it can be concluded in agree- 
ment with Nettleton [17] and Mulcahy and Smith [16] that the exact 
conditions, such as temperature, particle size, oxygen partial pres- 
sure, and ambient gas velocity, under which the chemical process or 
the mass transfer process control the carbon combustion rate, are un- 
known. 

0. Experimental Work 

1. Functional Dependencies - Experimental investigations to de- 
termine the effects of the surface temperature, oxygen partial pres- 
sure, and gas velocity and temperature on the combustion rate of car- 
bon have been conducted. Parker and Hottel [14] ran experiments with 
2.5 cm diameter rods and found that the specific combustion rate was 
proportional to the 0.37 power of the mass velocity of the gas. 

Kuchta, Kent, and Damon [15], using carbon cylinders in air with a 
velocity range of 8.6 to 165 m/sec and an air temperature range of 
900* to 1200* C, found that the absolute reaction rate was proportion- 
al to the 0.47 power of the gas velocity. Froberg [19], using 1.27 cm 
diameter spheres in air and oxygen flows of 0.1 to 1.0 cm/sec and 
above 760° C, found that the reaction rate was directly proportional 
to the 0.5 power of the gas velocity. Tu, Davis, and hottel [10], 
using 2.5 cm diameter spheres in the temperature range of 680* C to 
1430° C and a flow rate of 3.51 to 54.3 cm/sec, found the combustion 
rate varied as the 0.4 to 0.5 power of the mass velocity. Matsui, 
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Koyama, and Uehara [20] impinged a jet of oxygen on a flat carbon disk 
and found that the combustion rate was proportional to the 0.5 power 
of the velocity gradient for surface temperatures above 1227* C. All 
investigators agree that the combustion rate increases with increasing 
gas velocity when the process is controlled by mass transfer of oxi- 
dant to the surface. The combustion rate was found to increase pro- 
portionally to the mass velocity of the gas raised to the 0.37 to 0.7 
power. 

Investigators agree that when the combustion process is controlled 
by the chemical process, the combustion rate is strongly dependent on 
the surface temperature, and increases with increasing surface temper- 
ature. However, when mass transfer controls the reaction process, the 
combustion rate was found to be less dependent on temperature. 

Kuchta, Kant, and Damon [15] found that for a flow rate of approxi- 
mately 97.5 m/sec, the combustion rate increased proportionally to the 
0.20 to 0.25 power of the air temperature, which they state is of the 
order of magnitude set by theory. They also state that the effect of 
surface temperature is greater than that of the air temperature, but 
they were unable to quantitatively determine its effect. Tu, Davis, 
and hottel [10] found that the combustion rate varied with the arith- 
metic mean boundary layer temperature raised to the 0.6 to 1.1 power, 
which is within the limits of their theory. 

Arthur [21] and Tu, Davis, and hottel [10] found that the com- 
bustion rate was proportional to the partial pressure of oxygen raised 
to a power slightly less than unity, however, Rheed and Wheeler [22] 
found that the combustion rate was directly proportional to the par- 
tial pressure of oxygen. 
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Muleahy and Smith [*6] state that for mass transfer controlled 
combustion, a correction factor of approximately b was needed to bring 
the theoretical and experimental burning times of particles into 
agreement. When the chemical process was controlling the combustion, 
the prediction of the burning time was dependent on how accurate the 
pre-exponential factor, the activation energy, and the order of re- 
action were known. 

2. Primary Combustion Products - In the combustion of carbon, 
the primary products of combustion, i.e. , the products at the surface 
of the solid, are of major importance. The primary products can be 
CO, CO^, or a combination of CO and C0 2 . The products occurring 
at the surface of the particle can significantly affect the combustion 
rate, because one molecule of oxygen can remove either one carbon atom 

as C0 2 , or two carbon atoms as CO. lherefore, to accurately model 
the process and oetermine the reaction rate in carbon combustion, the 
primary products of combustion must be known. 

A summary of the work performed to determine the primary products 
of carbon combustion is presented in Table 3, and Table 4 presents the 
work done to determine the ratio CO/CO,, or C0/(C0 ♦ C0 2 ) formed 
during combustion. In all except four of the investigations, the gas 
analysis was performed on the end products of combustion. Experi- 
menters in some cases have used either a very fast flow rate to quench 
secondary reactions or a retardant to try to eliminate CO oxidizing to 
CO^. Four of the investigators used a probe to sample the gases in 
the combustion zone. Two of the investigators used uncooled probes 
and state that there may have been secondary reactions occurring in 
the probe. The other two investigators, Arthur [31] and Wicke and 


Table 3 - Combustion produ:*s 
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Wurzbacher [5], used a water cooled probe in order to quench any such 
secondary reactions. 

The experiments listed in Tables 3 and 4 were conducted under var- 
ious test conditions, which had an effect on the primary products 
formed. Arthur [21] found that the gas velocity had no effect on the 
CO/CO 2 ratios at low temperatures, but at high temperatures the 
CO/CO 2 ratio increased rapidly with gas velocity. The temperature 
at which this transition occurred was dependent on the type of car- 
bon. Phillips, Vastola, and Walker [28, 33] found that as the oxygen 
partial pressure increased, the C0/C0 2 ratio decreased. Lewis, 
Gilliland, and Paxton [27] agree with this result for partial pres- 
sures below 0.2 atm. Arthur [21] found the CO/CO 2 ratio to be inde- 
pendent of initial oxygen partial pressure in the range of 0.05 atm to 
0.25 atm. Phillips et al. [33] and Laine, Vastola, and Walker [29] 
found that the CO/CO 2 ratio increased with burn-off, but the effect 
was small except in the early stages where it increased substantial- 
ly. Lewis et al. [27] state that the C0/(C0 + CO^ ) ratio is rela- 
tively independent of the type of carbon while Walker, Rusinko, and 
Austin [34], Laine et al. [29], and Phillips et al. [33] state the 

CO/CO 2 ratio does depend on the nature of the carbon. Wicke and 
Wurzbacher [5] found that moisture in the air promoted the oxidation 

of CO to CO 2 . For moisture contents greater than 1%, CO hardly 
played a role in the oxidation process. 

Even though the experiments listed in Tables 3 and 4 were conduct- 
ed under various test conditions, it can be concluded that the 
CO/CO 2 ratio increases with temperature, which is generally accepted 
as being the case. And it can be seen that above 1100° C, a CO 2 
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reduction at the surface can occur and play a part in the combustion 
process. However, the temperature at which the primary products 

change from C0 2 to CO cannot be determined from these tables. 

Fendell [1], from his literature search, concludes that there are 
three zones of heterogeneou., combustion of carbon. A low temperature 
range (approximately 395* C to 700° C) in which C0 2 is the primary 
product of combustion; a mid-temperature range (approximately 700° C 
to 1100* C) in which CO and C0 2 are the products with a C0 2 re- 
duction at the surface occurring at approximately 895° C; and a high 
temperature range (approximately 1100° C to 2500* C) in which the 
following two different mechanisms for the combustion process have 
been proposed: (1) CO is formed at the carbon surface by direct oxi- 

dation and by the reduction of ambient CO^, and no gas phase re- 
action of CO occurs; and (2) C0 2 is reduced at the surface to form 
CO, which is burned in a gas phase reaction. 

3. Concentration and Temperature Profile Measurements - Wicke 
and Wurzbacher [5] experimentally obtained concentration profiles 
using a water cooled probe for vertical carbon tubes burning in dry 
and moist oxygen flows in the temperature range of 1100° C - 1230* C. 
Internally heated high purity spectal carbon tubes, 60 mm long, 15 mm 
outside diameter, and 3.5 mm wall thickness were used, bas samples 
were taken on the outside of the tube, 30 mm from either end. It was 
found that the concentration profiles obtained when dry oxygen (less 
than 2.6 x 10“^ percent water vapor present) was used were affected 
by the reduction reaction, C + C0 2 » 2C0 (Boudouard reaction), at 
the surface and the gas phase oxidation of CO to C0 2 . With a dry 
oxygen flow and at the lower surface temperatures, the 0 2 concentra- 


tion was nearly linear and extrapolated to near zero at the surface. 

As the surface temperature was Increased, the 0^ concentration pro- 
files became less steep and slightly curved. The CO profile was found 
to be independent of temperature and fell with distance due to homo- 
geneous burning of CO. The C0 2 concentration profile passed through 
a maximum, which moved to greater distances from the carbon surface 
and became more pronounced as the temperature increased. Also, in <...y 
oxygen, there was a 8 mm thick blue zone around the carbon tube, which 
was a result of CO being oxidized to C0 2 . With a 0.21% moisture 
content in the oxygen, the blue flame zone was reduced to a thickness 
of 1 or 2 mm. With a 1% moisture content in the oxygen, CO hardly 
played a role in the combustion process. 

Kisch [35], using spectrum line reversal and Ft Rh - 10% Pt 
thermocouples, was able to obtain temperature profiles for burning 
graphite. Graphite blocks, 20 x 20 x 50 mm, were internally heated 
and burned in dry oxygen flowing at 30 cm/ sec. Both the line reversal 
and thermocouple measurements showed that there was a maximum temper- 
ature between 1 and 2 mm in front of the burning graphite surface. 

For a 1127° C graphite surface temperature, the maximum gas phase tem- 
perature was 1627° C. 

Temperature profiles in front of small carbon rods, positioned 
vertically in a horizontal 6 cm diameter electric pipe furnace, were 
obtained by deGraaf [36] using the method of spectrum line reversal. 
The burning rates of the rods were changed by inclining the pipe furn- 
ace slightly. At a furnace temperature of 700° C, the carbon flared 
up and a blue flame zone formed around it. At 800° C, the blue flame 
zone was approximately 2.5 mm thick. A spectrum of the blue reaction 
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zone was recorded, and it was identical to the spectrum of a carbon 
monoxide-air flame. For a furnace temperature of 1000* C, it was 
found that as the burning rate increased, both the rod surface temper- 
ature and the difference between the maximum gas phase temperature and 
the rod surface temperature increased. The maximum temperature 
occurred approximately 1.25 mm from the surface. Gas phase temper- 
ature profiles were also obtained for different furnace temperatures. 
The difference between the maximum gas phase temperature and the rod 
surface temperature was found to decrease as the furnace temperature 
increased from 950* C to 1100* C. When the furnace temperature was 
increased from 1100* C to 1180* C, the difference between the maximum 
gas phase temperature and the carbon surface temperature increased. 

In this series of tests, the maximum gas phase temperature occurred 
between 1.0 mm to 1.5 mm from the carbon surface. Also, at the high- 
est furnace temperatures the maximum temperature differences were less 
pronounced than at the lower furnace temperatures. 

In summary, investigators have studied carbon combustion both 
theoretically and experimentally. The theoretical work has been con- 
cerned with the mathematical modeling of the combustion process, while 
the experimental work has been concerned with determining whether the 
combustion process is controlled by the mass transfer or by the chem- 
ical process. In addition, the effects of particle diameter, surface 
temperature, oxygen partial pressure, and gas velocity and temperature 
on the carbon combustion rate have been studied, and experiments have 
been conducted to determine the species concentration and temperature 
profiles. 


The wide variation in the results of the previously mentioned in- 
vestigations, along with their generally macroscopic nature, served as 
the motivating reasons for conducting the present study. It was felt 


- . .. . 


that a series of carefully conducted carbon combustion experiments 
under both zero and normal gravity conditions would more clearly quan- 
tify convective effects, illustrate a possible mechanism dependence on 
surface position, and provide data more in line with the assumptions 
made in the mathematical models. 


III. ZERO GRAVITY EXPERIMENTATION 


A. lest Facility 

The experimental data for this part of the study were obtained in 
the Zero-Gravity Facility at the Lewis Research Center. A schematic 
diagram of this facility is shown in Figure 5. The facility consists 
of a cone rete-1 'red 8.5 meter (28 ft) diameter shaft that extends 155 
meters (510 ft) below ground level. A steel vacuum chamber, 6.1 
meters (20 ft) in diameter and 143 meters (470 ft) high, is contained 
within the concrete shaft. The pressure in this vacuum chamber is 
reduced to 13.3 newtons per square meter (1.3 x 10“^ atm) by util- 
izing the center's wind tunnel exhaust system and an exhauster system 
located in the facility. 

The ground-level service building has, as its major elements, a 
shop area, a control room, and a clean room. Assembly, servicing, and 
balancing of the experiment vehicle are accomplished in the shop 
area. Tests are conducted from the control room, Figure 6, which con- 
tains the exhauster control system, the experiment vehicle predrop 
checkout and control system, and the data retrieval system. 

1. Mode of Operation - In the Zero-Gravity Facility, the experi- 
ment vehicle is allowed to fall freely from the top of the vacuum 
chamber, resulting in nominally 5 seconds of free-fall time. By 
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Figure b. - Schematic diagram of zero-gravity facility. 
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free-fall is meant that there are no guide wires* electrical lines* 
etc., connected to the vehicle. Therefore, the only force (aside from 
gravity) acting on the freely falling experiment is due to residual 
air drag. This results in an equivalent gravitational acceleration 
acting on the experiment which is estimated to be on the order of 
10“^ g maximum. 

2. Recovery System - After the experiment vehicle has traversed 
the total length of the vacuum chamber, it is decelerated in a 3.6 
meter (12 ft) diameter, 6.1 meter (20 ft) deep container which is lo- 
cated on the vertical axis of the chamber and filled with small pel- 
lets of expanded polystyrene. The deceleration rate (averaging 32 
g's) is controlled by the flow of pellets through the area between the 
experiment vehicle and the wall of the deceleration container. The 
deceleration container is shown in the photograph of Figure 7. 

B. Experimental Vehicle 

The experiment vehicle used to obtain the data for this study is 
shown in Figure 8. The overall height (exclusive of the support 
shaft) was 3.0 meters (9.85 ft), and the largest diameter was 1.06 
meters (3.5 ft). The vehicle consisted of a telemetry section con- 
tained in the aft fairing and an experimental section housed in the 
cylindrical midsection. 

1. Telemetry System - The on-board telemetry system which was 
used to record pressure data was a standard Inter-Range Instrumenta- 
tion Group (IRIG) FM/FM 2200-megahertz telemeter. It was used during 
a test drop to record up to 18 channels of continuous data. The sys- 
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tem had a frequency range up to 2100 hertz. The telemetered data were 
recorded on two high-response recording oscillographs located in the 
control room. 

2. Experimental Section - The details of the experimental sec- 
tion are shown in Figure 9. The spherical tanks at the top section 
were part of a gaseous nitroqen fire extinguishment system which was 
used as a safety precaution at the completion of each test. The bulk 
of the experimental equipment was contained in the central portion of 
the drop module and consisted of a combustion chamber with a carbon 
rod holder, two high-speed cameras, power supplies, electrical control 
box, and plumbing used to fill and evacuate the chamber. 

A cross-sectional schematic of the combustion chamber is shown in 
Figure 10. The combustion chamber was a 78.74 cm (31 in.) long by 
39.62 cm (15.6 in.) diameter cylinder with a 19.81 cm (7.8 in.) radius 
spherical cap. The internal volume of the chamber was approximately 
1*1 x 10 5 cm^ (6,900 in.^). The walls of the chamber were 
0.478 cm (0.188 in.) t! . k .tainless steel. There were four ports for 
electrical and instrument’. icn wires. Three of the ports were 12. 7u 
cm (5 in.) from tne bottom of the chamber and were 90° apart. The 
fourth port was 66.04 cm (26 in.) from the bottom of the chamber. 

There were also two viewing ports that were used to photograph the 
combustion process. One viewing port was 12.70 cm (5 in.) from the 
bottom of the chamber and the other one was at the top of the spher- 
ical cap. The chamber was split 25.40 cm (10 in.) from the bottom so 
that the top part of the chamber could be removed for replacement of 
the experiment. The drop module was equipped with rails and rollers 
so the top portion of the chamber could be easily moved. 
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The carbon rod burning apparatus was inside the combustion chamber 
as shown in Figure 11. The carbon rod holder. Figure 12, consisted of 
a 22.86 cm (9 in.) square by 0.478 cm (0.188 in.) thick stainless 
steel plate with two 0.478 cm (0.188 in.) diameter bolts 3.81 cm (1.5 
in.) long in it. The two bolts were 7.62 cm (3 in.) apart. Lead 
wires were held to the bolts by tie downs. At the end of the lead 
wires were junctions that connected to tungsten electrodes. These 
electrodes were pushed into the ends of the carbon rods. The diameter 
of the electrodes varied depending upon the diameter of the carbon rod 
used in the test. When a 0.615 cm (0.242 in.) diameter or 0.457 cm 
(0.180 in.) diameter carbon rod was used, the electrodes were 0.160 cm 
(0.063 in.) in diameter. When a 0.305 cm (0.120 in.) diameter carbon 
rod was used, the electrodes were 0.102 cm (0.040 in.) in diameter. 

Two cameras were used to photograph the combustion process. One 
camera was located at the side viewing port and photographed the end 
of the carbon rod. This camera ran at 400 frames/sec. The other 
camera was positioned at the top of the combustion chamber and photo- 
graphed the length of the rod from the top. This camera ran at 200 
frames/ sec. A time mark every 0.01 sec was put on the edge of the 
film. Also, a continuous mark was put on the edge of the film before 
the package was released. At release, this mark was cut off giving a 
reference for the start of zero-gravity. 

Power to ignite the carbon rod was supplied by rechargeable lead 
acid battery packs. For the tests using a 0.615 cm (0.242 in.) diam- 
eter rod, 48 volts were used and the current draw was approximately 
220 amps. For tests using a 0.457 cm (0.180 in.) diameter rod, 32 
volts were used with a currant draw of approximately 140 amps. Fin- 





43 


ally, for tests using a 0.305 cm (120 in.) diameter rod, 16 volts were 
used with a current draw of approximately 50 amps. A separate power 
supply using rechargeable nickel-cadmium batteries was used to power 
the other equipment in the drop vehicle. 

C. Test Materials 

The carbon rods used in the tests were Ultra-Carbon grade U-5 pe- 
troleum coke-based spectrographic carbon rods and were, as stated, 
0.615 cm, 0.457 cm, or 0.305 cm in diameter. The rods were cut into 
5.08 cm (2.0 in.) lengths. Holes of the diameter of the electrodes to 
be used for a certain sized rod were drilled 0.396 cm (0.156 in.) into 
the ends of the rod. The atmosphere for the tests was 99.996% pure 
oxygen with a maximum water vapor content of 0.5 ppm. 

D. Test Procedure 

Prior to each run, the combustion chamber was wiped clean. Before 
the carbon rods were inserted into the holder, the holes in the ends 
of the rods were filled with silver brazing powder. The powder was 
used to reduce the contact resistance between the rod and the elec- 
trode in order to provide a more uniform heating of the rod. After 
the rod was on the electrode, it was lined up along the axis of the 
camera that photographed its end. The combustion chamber was then 
closed and sealed. The chamber was evacuated and purged with dry ni- 
trogen three times before it was filled with dry oxygen. 


The proper voltage for ignition of the size rod used in the test 
was set, along with the length of time power would be supplied to the 
electrodes. Power to ignite the carbon rods was set to go on 0.045 
sec after the vehicle was released. This prevented free convection 
currents due to the heated rod from being set up inside tho combustion 
chamber. 

The vehicle was positioned at the top of the vacuum chamber as 
shown in Figure 13. It was suspended by the support shaft on a 
hinged-plate release mechanism. During the vacuum chamber pump-down 
and prior to release, monitoring of experiment vehicle systems was 
accomplished through an umbilical cord attached to the top of the sup- 
port shaft. Electrical power was supplied from ground equipment, with 
the system being switched to internal power a few minutes before re- 
lease. The vehicle was released by pneumatically shearing a bolt that 
was holding the hinged plate in the closed position. No measurable 
disturbances were imparted to the experiment vehicle by this release 
procedure. 

The total free fail test time obtained was 5.16 sec. During the 
test drop, the vehicle's trajectory and deceleration were monitored on 
closed circuit television. Following the drop, the vacuum chamber was 
vented to the atmosphere and the experiment was returned to ground 
level . 

E. Data Reduction 

The burning rods were photographed on high-speed film which was 
examined on a motion picture analyzer. Flame diameters were measured 
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as a function of time from the film showing the end viewof the rod. 

No data on the flame could be obtained from the film showing the 
length of the rod from above because the brightness of the glowing rod 
obscured the pale blue flame that surrounded the rod. Also, the dis- 
tance the camera was away from the rod made photographing of the flame 
difficult even though this camera ran half as fast as the camera 
photographing the end of the rod. 

F. Experimental Results and Discussion 

Tests burning 0.615 cm (0.242 in.), 0.457 cm (0.180 in.), and 
0.305 cm (0.120 in.) diameter by 5.08 cm (2.0 in.) long carbon rods in 
3.45 x 10^ N/m 2 (5 psia), 6.89 x 10** N/m 2 (10 psia), 1.03 x 
10 5 N/m 2 (15 psia), and 1.38 x 10^ N/m 2 (20 psia) dry oxygen 
environments were conducted under zero-gravity conditions. The rods 
were ignited by passing an electric current through them for a set 
time. The average times for the 0.615 cm, 0.457 cm, and 0.305 cm rods 
were 1.04 sec, 0.823 sec, and 1.348 sec, respectively. The rods 
were allowed to burn, supported by their own combustion, for the re- 
mainder of the zero-gravity time. In some tests, the power to the 
carbon rod was not left on long enough to completely ignite the rod. 

In all except one of the incomplete ignition tests, there was no flame 
but the rod glowed and decreased in intensity throughout the duration 
of the test, indicating that the reaction was being quenched. In the 
other incomplete ignition test which involved a 0.615 cm (0.242 in.) 
diameter rod in 1.38 x 10^ N/m 2 (20 psia) oxygen, the flame 
flickered and went out as the rod cooled during the test. In all 
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other tests where sustained burning occurred in zero-gravity, a blue 
flame surrounded the rod. Photographs of the end of the burning rod 
'or two tests are shown in Figures 14 and 15. In Figure 14, a 0.615 
cm (0.242 in.) diameter rod in 1.38 x 10^ N/m 2 (20 psia) oxygen 
is shown, and the same diameter rod in 3.45 x 10 4 N/m 2 (5 psia) 
oxygen is shown in Figure 15. From these photographs it can be seen 
that a symmetrical blue flame extending from the surface surrounded 
the rod, indicating that the CO formed at the surface was being burned 
to CO 2 in a gas phase reaction throughout the region. This means 
the carbon rod was burning according to one cf the high temperature 
models that allows a gas phase reaction. Because the flame extended 
from the su. f ace of the rod, the carbon was burning according to the 
models presented by Hugo, Wicke, and Wurzbacher [4] or Caram and 
Amundson [6], If the flame did not extend from the surface, that is, 
if it stood off from the surface, the carbon would be burning accord- 
ing to the models presented by Spalding [7] or van der Held [8]. To 
obtain the surface temperature, an optical technique would have had to 
ha used, but this was not possible in the zero-gravity facility. 

1. Measured Data - The data obtained for the 0.615 cm (0.242 
in.), 0.457 cm (0.180 in.), and 0.305 cm (0.120 in.) diameter rods, at 
four different pressures of 1.38 x 10 5 N/m 2 (20 psia), 1.03 x 
10 5 N/m 2 (15 psia), 6.89 x 10 4 N/m 2 (10 psia), and 3.45 x 
10^ N/m 2 (5 psia) are presented in Figures 16, 17, and 18, re- 
spectively. The data presented in the figures were averaged for the 
tests that were conducted more than once. These figures show the 
effect of pressure on the ratio of flame diameter to rod diameter. For 
the three rod diameters, the ratio of flame diameter to rod diameter 


9 


fiour? 14 CJfiun rod <0. 6IS cm dum. i burning in 
1. W*.(r N/nr (?Cjhiaionygen environment. 
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Figure IV Carbon rod IU 61* cm dum. » burning i 
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Figure 18. - Effect of pressure on the ratio of flame diameter to rod diameter for a 3® 


increased as the pressure decreased. Also, the ratio of flame diam- 
eter to rod diameter increased less rapidly as the burning time in- 
creased, From Figures 16 and 17, it can be seen that for the 0.61b cm 
(0.242 in.) and 0.457 cm (0.100 in.) diameter rods at times greater 
than 4 sec, the data generally tends to level off which seems to in- 
dicate the ratio may have been approaching a steady value. 

The average data for the 1.38 x 10^ N/m 2 (20 psia), 1.03 x 
10 5 N/m 2 (15 psia), 6.89 x 10 4 N/m 2 (10 psia), and 3.45 x 
10 4 N/m 2 (5 psia) runs for the three rods used are presented in 
Figures 19, 20, 21, and 22. From these figures it can be seen that 
the ratio of flame diameter to rod diameter decreased as the rod diam- 
eter increased. Figures 19 and 20 show that the ratio of flame diam- 
eter to rod diameter tends to level off for times greater than 4 sec 
for the 0.615 cm (0.242 in.) and 0.457 cm (0.180 in.) diameter rods 
for pressures of 1.38 x 10 5 N/m 2 (20 psia) and 1.03 x 10^ 

N/m 2 (15 psia) . 

Figures 16 through 22 show that the ratio of flame diameter to rod 
diameter increased less rapidly as the burning time increased. As the 
flame diameter increases, the area for diffusion as well as the volume 
over which the flame occurs increases, making more oxygen molecules 
available for the gas phase reaction. This means the flame will grow 
less rapidly with time because there is more oxygen available for the 
reaction. 

2. Statistical Analysis - The data obtained from each zero-grav- 
ity test was curve-fitted numerically using an exponential fit. The 
equations obtained along with the sum of residuals squared (£ resid- 
uals 2 ), the total sum of squares (total SSQ), R squared (R 2 ) val- 
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Figure W. - Effect of rod diameter on the ratio of flame diameter to rod diameter for 1.38x10? N/m 2 (20 psia) pressure. 
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ues, the number of data points (n) used in the curve fit, the variance 
around the fitted line (MSE) and the standard deviations (S) are pre- 
sented in Table 5. The sum of residuals squared is 2(d* _ d*) 4 , 

A 

where d* is the value obtained from the curve fit and d* is the actual 
ratio of flame diameter to rod diameter. The total sum of squares is 
-(dT _ d*) 2 and is equal to 2(d? - d?) 2 + £(<■}? - d*) 2 
where d* is the mean value of d*. The R 2 is equal to (total 
SSQ - 2 residuals 2 )/total SSQ or - d*) 2 / £(d* - d*) 2 . 

According to Reference [37], the R 2 value is a measure of the pro- 
portion of the total variation about the mean explained by regres- 
sion. Therefore the R value is a measure of the goodness of fit; 
if R = 1 there is a good fit and if it is zero there is no fit. 

The R values obtained for the curve fits ranged from 0.974 to 0.577 
which means, in general, that the fits are fairly good. The variance 
around the fitted line is the sum of the residuals squared, divided by 
the number of points used in the fit. The standard deviation is the 
square root of the variance around the fitted line. 

The experimental data along with the curve fits are presented in 
Figures 23 to 34. These figures show that the curve fits do fit the 
data well. Also, Figures 26, 27, 28, 29, 31, 32, and 33 show that, in 
general, the variation within a test was less than the variation be- 
tween tests that were conducted with the same diameter rods and at the 
same pressure. The differences between tests are attributed to the 
surface temperature of the rod, which could not be controlled or meas- 
ured. If sufficient time were available so that the rods could have 
reached an equilibrium burning temperature, they would be expected to 
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Figure 33. - Comparison of experimental data and curve fit for a 615 cm (a 242 in . ) diameter rod burned in L 03x11? film 2 (15 
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have the same ratio of flame diameter to rod diameter for tests con- 
ducted at the same pressure and with the same diameter rods* 

The coefficients from the curve fit equations were used to obtain 
a correlation equation for all the data. The equation obtained is: 

3 * • o(l - e“^l t- ^2) 

where a = 1.943 - 6.599(d - 0.18) - 0.03517 (P - 12.5) ♦ 

46 . 69 ( d - 0.18)2 
■ 0.687 

t> 2 * 1-03 5 + 4.335(d - 0.18) * 0.0278(P - 12.5) - 
119. 4{d - 0.18)2 


The units of d (rod diameter), P (pressure), and t (time) are in., 
psia, and sec, respectively. 

The standard deviations of a, and P 2 about the regression 
line are 0.201, 0.219, and 0.196, respectively. The coefficients in a 

and ^re those that had a T ratio greater than 1.98. Only the 
46.69 coefficient in o had a value this low, with all others being 
greater than two. The T ratio, which is the coefficient divided by 
the standard deviation, indicates if there is a relationship between 
the vu^iables. The larger the value of T, the better the relation- 
ship. If T is les. than 2, there is no relationship, and if it is 
greater than 2, there is a relationship. Because only the first coef- 
ficient in had a T value greater than 2, the mean value of 
was used. 
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The data along with the correlation equation are presented in Fig- 
ures 35 through 38. The data on each of the figures is presented with 
a constant pressure and variable rod diameter. In Figure 35, it can 
be seen that at 3.45 x 10^ N/m^ (5 psia), the correlation equation 
fits the 0.305 cm (0.120 in.) diameter rod data, but is slightly 
higher than the 0.457 cm (0.180 in.) and the 0.615 cm (0.242 in.) di- 
ameter rod data. Figure 36 shows that at 6.89 x 10^ HI m 2 (10 
psia), the correlation equation fits the 0,305 cm (0.120 in.) and the 
0.457 cm (0.180 in.) diameter rod data, while it is slightly higher 
than the 0.615 cm (0.242 in.) diameter rod data. In Figure 37, it can 
be seen that at 1.03 x 10^ N/m^ (15 psia), that the 0.305 cm 
(0.120 in.) diameter rod data does not exactly follow the correlation 
equation, and the correlation equation is slightly lower than the 
0.457 cm (0.180 in.) diameter rod data. The correlation does fit the 
0.615 cm (0.242 in.) diameter rod data for times grater than 2.5 sec- 
onds. Figure 38 shows that at 1.38 x 10 5 N/m 2 (20 psia), the cor- 
relation equation fits the 0.615 cm (0.242 in.) and the 0.457 cm 
(0.180 in.) diameter rod data, but the 0.305 cm (0.120 in.) diameter 
rod data falls below the correlation equation. From Figures 35 
through 38, it can be concluded that, in general, the correlation 
equation predicts the ratio of flame diameter to rod diameter (d*) 
within approximately +0.1, except for the 1.38 x 10 5 N/m 2 (20 
psia) data for the 0.305 cm (0.120 in.) rod where the correlation 
equation and the ratio of flame diameter to rod diameter vary by as 
much as +0.25. 








•or rods burned in L C8*1C? (15 psia) oxygen. 
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6. Conclusions for Zero Gravity Tests 

A blue flame surrounding the carbon rod under zero-gravity condi- 
tions showed that a gas phase reaction occurred, in which carbon mon- 
oxide was oxidized to carbon dioxide. This means that the carbon rod 
was burning according to one of the high temperature models that 
allows a gas phase reaction. Because the flame extended from the sur- 
face of the rod, the carbon would seem to be burning according to the 
models presented by Hugo, Wieke, and Wurzbacher [4] or Caram and 
Amundson [6], If the flame had not extended from the surface, that 
is, if it stood off from the surface, the carbon would be burning 
according to the models presented by Spalding [7] or van der held 
[8], It was not possible to obtain the surface temperature in the 
zero-gravity tests. From the zero- 3 ravity data, it was found that the 
ratio of flame diameter to rod diameter increased as the pressure de- 
creased, and decreased as the rod diameter increased. The data was 
used to obtain the correlation equation: 

d* = a( 1 - e^l 1 -^) 

where a = 1.943 - 6.599(d - 0.18) - 0.03517 (P - 12.5) + 
46.69(d - 0.18) 2 

and pj , 0.687 

t> 2 « 1.015 ♦ 4.335(d - 0.18) ♦ 0.0278(P - 12.5) - 
119. 4(d - 0.18) 2 


The units of d (rod diameter), P (pressure), and t (time) are in., 
psia, and sec, respectively. In general, this equation predicts the 
ratio of flame diameter to rod diameter (d*) within approximately 
+0.1, except at 1.38 x 10^ N/m^ (20 psia) for a rod diameter of 
0.305 cm (0.120 in.), where the deviation approaches +0.25. The in- 
formation obta.ned from this eAperiment can be used to properly design 
an in-space experiment, such as for Spacelab, where long-term 
zero-gravity information on the combustion process can be obtained. 

The long-term zero-gravity data should include surface temperatures 
and temperature and concentration profiles, which are needed to accu- 
rately model the zero-gravity combustion process. The concentration 
profiles would confirm which carbon-burning model is correct. If : ,he 

O 2 concentration extended to the surface of the rod, the models pre- 
sented by Hugo, Wicke, and Wurzbacher [4] or Caram and Amundson [6] 

would be correct. If the O 2 concentration goes to zero before the 
surface of the rod is reached, the models presented by Spalding [7] or 
van der held [8] would be correct. 


IV. NORMAL GRAVITY EXPERIMENTATION 


A. Apparatus 

The experimental system used to study the products of combustion 
from a burning carbon rod in normal gravity consisted of a carbon rod 
igniter and a modulated beam mass spectrometric sampler. 

The carbon rod igniter and holder is shown in Figure 39. In all 
of the normal gravity tests, 0.615 cm (0.242 in.) diameter by 5.08 cm 
(2.0 in.) long carbon rods were used. The carbon rod was held by two 
tungsten electrodes, which were 0.160 cm (0.063 in.) in diameter and 
extended through two copper poles. The center of the carbon rod was 
6.1 cm (2.4 in.) above the flow straightener. A glass cylinder with a 
10.80 cm (4.25 in.) inside diameter and a height of 11.43 cm (4.5 in.) 
covered the carbon rod holder. 

The rods were burned at atmospheric pressure in oxygen with a 

water vapor content of 6 ppm. A volumetric flow rate of 0.15 1/sec 
3 

(0.005 ft /sec) was set by a needle valve and measured with a cali- 
brated rotameter. This volumetric flow rate, which had to be used to 
make up for the gases drawn off by the sampling probe, corresponded to 
a 4.35 *'m/sec (0.159 ft/sec) flow rate past the rod. The oxygen 
entered through the bottom of the lower section of the carbon rod 
holder, which was a 7.62 cm (3.0 in.) long cylinder with a 6.v% cm 
(2.7 in.) inside diameter. The oxygen flowed through a screen and 
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then through a section containing plastic balls, 0.572 cm (0.225 in.) 
in diameter, in order to spread the flow. Finally, the oxygen passed 
through a Hastelloy X honeycomb flow straightener. The flow straight- 
ener had honeycomb openings 0.08 cm (0.031 in.) on a side with 135 
openings/cm^ (871 openings/ in. ^), giving an open area of 29.31 
cm^ (4.543 in.^). The burner was support'd by mechanical devices 
which facilitated micrometer movement in three mutually perpendicular 
directions. Power to ignite the carbon rod was supplied by a 48 volt 
power supply. 

Direct mass spectrometric analysis of the species in the atmo- 
snheric pressure flame surrounding the carbon rod was accomplished 
with a modulated molecular beam mass spectrometer sampler, shown in 
Figure 40. The sampler is described in more detail in Reference [38] 
and is shown schematically in Figure 41. Atmospheric pressure sam- 
pling was accomplished by drawing a gas sample through either a 0.20 
mm (0.008 in.) diameter Pt-10%Rh orifice shown in Figure 41, or 
through a quartz sampling probe with a 0.20 mm (0.008 in.) diameter 
orifice shown in Figure 42. Quartz sampling probes with orifice diam- 
eters of 0.03 mm (0.001 in.), 0.13 mm (0.005 in.), and 0.25 mm (0.010 

in.) were also tried. The smaller probes (0.03 mm and 0.13 mm) plug- 

ged up with small flakes from the carbon surface and the 0.25 mm ori- 
fice was so large that the sampling system could not be pumped down. 
The Pt-10% Rh orifice was used to obtain samples at the top of the rod 
(the 0° position), and the quartz probe was used to obtain samples 
at positions of 45° and 90° from the top of the rod. 

To sample the center portion of the flow from the sampling ori- 
fice or probe, the stream was passed through a skinner cone with an 
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Figure 41. - Schematic of spectrometric sampler. 
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orifice diameter of 0.81 mm (0.032 in.) and an included angle of 
60°. The distance between the sampling orifice or the top plate of 
the quartz sampling probe and the skimmer orifice was 3.18 cm (1.25 
in.). For this physical configuration with nominal pumping speeds 
indicated in Figure 41, first and second stage pressures were approxi- 
mately 1.5 x 10”3 torr and 8 x 10"^ torr, respectively, when sam- 
pling room temperature, atmospheric pressure gas. When the gas temp- 
erature was increased to 1000* C, these pressures were 7 x 10“** torr 
£ 

and 1 x 10 torr, respectively. Third and fourth stage pressures 

7 ft 

were always iess than 10 and 10 torr, respectively, for all 
sampling conditions. Stage I pressures were read with a capacitance 
manometer, and the pressures in the other three stages were read with 
ion gauges. The molecular beam from the skimmer was chopped by a 
motor driven two-toothed chopper wheel located in Stage II. A chop- 
ping frequency of 150 Hertz was used and a reference signal at this 
frequency was derived from a light bulb and photodiode coupled to the 
chopper wheel. The chopped molecular beam then passed to the electron 
bombardment ion source of a quadrupole mass spectrometer. The quad- 
rupole filter with 1.6 cm diameter poles had a mass range extending to 
over 600 AMU. A Channeltron electron multiplier was employed to mul- 
tiply the ion current output of the ion seurce-quadrupole filter. Two 
channels of current output were measured as a function of quadrupole 
filter tuning. One channel measured the total chopped ion current, 
and the other channel measured only the component of the ion current 
signal in phase with the chopper. The second channel was driven by a 
lock-in amplifier-phase sensitive detector system tuned to the chopper 
frequency reference signal. The gaseous species in the flame were 
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measured by recording the in-phase component of the ion current for 
respective values of the mass to charge ratio. The accuracy of the 
mole fractions obtained is uncertain because the system was not cali- 
brated for such mass dependent factors as quadrupole mass filter 
transmission, multiplier gain, or Mach number focusing. However, be- 
cause the mass range used was narrow, 26 to 46 AMU, the errors associ- 
ated with mass dependent factors are estimated to be less than 10 per- 
cent. 

B. Test Procedure 

The sampling orifice to be used in the test was attached to the 
mass spectrometric sampler. The sampling apparatus was pumped down 
and put into operation and the sampling orifice position was adjust- 
ed. The adjustment was made while ambient laboratory air was being 
sampled. The mass filter was set on the maximum of the nitrogen peak 
and the electrometer output on the recorder was observed while adjust- 
ing the sampling orifice centering screws. Adjustments were made to 
maximize the nitrogen signal. The quadrupole mass filter range was 
set for the range of 26 to 46 AMU with a rate of 0.2 sec/AMU, so 
spectra were recorded every four seconds. 

The noles in the ends of the carbon rods were filled with silver 
brazing powder before the rod was put into the carbon rod holder in 
order to reduce the contact resistance between the electrodes and the 
rod. 

For the tests where gas samples were taken at the top of the rod, 
the glass cylinder was moved up to the plate holding the orifice, and 
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the carbon rod was positioned 0.69 cm (0.27 in.) below the sampling 
orifice. After the rod was ignited, it was moved toward the sampling 
orifice. At leas . two spectra were taken at each position to which 
the carbon rod was moved. The distance the carbon roo was moved was 
measured with a dial indicator that had increments of 0.003 cm (0.001 
in.). For the tests where the quartz probe was used to sample the 
gases, the probe was set 0.025 cm (0.010 in.) from the surface and 
left at this position while the rod burned. Pecause the quartz sampl- 
ing probe was lower than the Pt-10% Rh probe, the glass cylinder was 
open at the top. To eliminate ambient air from entering the carbon 
rod holder, aluminum foil with a slot in it for the sampling probe was 
placed on top of the glass cylinder. 

Before a rod was ignited, the carbon holder was purged with oxy- 
gen. The gas species in the carbon holder assembly were monitored 
with the mass spectrometer. The rod was ignited when oxygen was the 
only gas present in the assembly, with ignition obtained by passing an 
electric current through the rod. Power was supplied by a 48 volt 
power supply, and the power was left on until a flame appeared around 
the carbon rod. 

C. Data Reduction 

The ion intensities were corrected for relative ionization cross 
sections. Cross sections for atoms were taken from Mann [39] and est- 
imated for molecules as 0.75 times the sum of the atomic values, be- 
cause it was impossible to photograph the center of the rod where the 
samples were taken while the carbon holder assembly was in position in 
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the mass spectrometer system, the carbon holder assembly was pulled 
out and three rods were burned in it and photographed from above and 
from the side. One of the rods had a flame only over the top sur- 
face. Of the other two rods, only one was photographed until it com- 
pletely burned. A side view of the carbon rod burning is shown in 
Figure 43. The respective diameters of this rod as a function of time 
for the two viewing positions are presented in Figure 44. It can be 
seen that the sides of the rod (top view) burned faster than the top 
and bottom combined (side view). Also, the side view change of posi- 
tion of the horizontal centerline of the rod from its original posi- 
tion is shown in Figure 45. This figure shows that the bottom of the 
rod burned faster than the top. Based on the data in Figures 44 and 
45, the rod became elliptical as it burned. To determine the probe 
distance from the surface, the diameters of the rod at the time the 
spectra were taken were used to determine the surface position of the 
ellipse. The change in surface position between the original carbon 
surface and the elliptical surface was added to the original probe 
position to get the horizontal distance between the probe and the sur- 
face. 


D. Experimental Results and Oiscussion 

Concentration profiles of the combustion products of a carbon rod 
as a function of distance from the surface were obtained for initial 
sampling probe positions of 90°, 45°, and 0°, with 0° being 
the top of the rod. As stated previously, the lower part of the rod 
burned faster than the top of the rod, so that the angle of the probe 








initially positioned at 45° and 90° changed as the rod burned. In 
Figure 46, the concentrations of 0 2 , CO?, and CO as a function of 
distance from the surface for the 90° quartz sampling probe are pre- 
sented. The 0 2 concentration increases with distance from the sur- 
face, while the CO concentration decreases with distance from the sur- 
face. The C0 2 concentration at first increases slightly and then 
decreases with distance. The concentration of CO is lower than the 
concentration of C0 2 except near the surface, that is, less than 
0.05 cm (0.02 in.) from the surface. Also, at this distance the con- 
centration of C0 2 can be seen to level off. None of the concentra- 
tions seem to go to zero at the surface, which indicates that the sur- 
face conversion process at this position was chemical process or, ne- 
glecting the chemisorption and desorption steps, chemical reaction 
rate controlled rather than mass transfer controlled. 

The concentration profiles obtained with the quartz probe at the 
45° position are presented in Figure 47. Again, the 0 2 concentra- 
tion increases with distance from the surface, while both the C0 2 
and CO concentrations decrease with distance. The 0 2 concentration 
is lower than the C0 2 and CO concentrations very close to the sur- 
face, and the CO concentration is always lower than the C0 2 concen- 
tration. The concentrations, if extrapolated to the surface, do not 
appear to equal zero, but the 0 2 profile does indicate that it would 
have a low value. This again indicates a chemical reaction rate con- 
trolled surface conversion at the 45° probe position. 

The concentration profiles obtained with the Pt-10% Rh probe at 
the top of the rod are presented in Figure 48. At this position, as 
at the other positions, the 0 2 concentration increases with distance 
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while the CO concentration decreases. The C0 2 at first increases 
with distance from the surface, reaching a maximum at approximately 
0.36 cm (0.14 in.) from the surface, and then starts to decrease. The 
concentration of C0 2 is higher than that of the other species except 
for 0 2 at large distances from the surface, greater than 0.64 cr.i 
(0.25 in.). The rate of CO concentration decrease is more rapid near 
the surface than away from the surface. The concentration of 0 2 is 
lower than the concentration of CO or C0 2 near the surface, approxi- 
mately 0.20 cm (0.08 in.) from the surface. Also, at 0.64 cm (0.25 
in.), there is a rapid increase in the 0 2 concentration and a cor- 
responding rapid decrease in the C0 2 concentration, indicating that 
the ambient oxygen in the carbon rod holder assembly is mixing with 
the combustion products. Data at distances less than 0.1 cm from the 
surface for the 0° probe position could not be obtained, because by 
the time sampling this close to the surface could be done, the rod 
diameter was so small that ambient gas was being drawn into the probe. 

The 0 2 concentration as a function of distance from the surface 
for the three sampling positions is presented in Figure 49. The 0 2 
concentration decreases with position from 90° to 0°. The O 2 
profile tends to reach the same value for the 90° and 45° probe 
positions at an approximate distance of 0.23 cm (0.09 in.) from the 
surface. If the concentrations are extrapolated to the surface, it 
can be seen that the surface concentration decreases considerably from 
the 90° position to the 45° position. It would be difficult to 
extrapolate the 0° position to the surface, but if the trend set up 
by the 90° and 45° probe positions holds, it is very likely that 
at the 0° position the surface concentration would go to zero. This 
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DISTANCE FROM SURFACE, cm 

Figure 49. - O? concentration as a function of distance from surface for various gas sampling probe 
positions. 


indicates that the oxygen conversion process at the surface is chem- 
ical reaction rate controlled at both the 90° and 45° probe posit- 
ions, and that it is very possible that the process at the 0° probe 
position is mass transfer controlled. 

The CO 2 concentration as a function of distance from the sur- 
face for the three probe positions is presented in Figure 50. The 

CO 2 concentration is higher ior the 45° probe position than it is 
for the 90° probe position. The concentrations for these two posi- 
tions tend to become equal to each other away from the surface, that 
is, at approximately 0.23 cm (0.09 in.) from the surface. The concen 
trations at the 0° probe position are greater than those at the 
other two probe positions. The positive slope of the CO 2 concentra- 
tion profiles for the 90° and 0° probe positions indicates that 
the CO 2 is reacting at the surface to produce CO. The increase of 
C0 2 with position from the 90° to 0° indicates that the combust- 
ion products are being convected upward. 

The CO concentration as a function of distance from the surface 
for the three sampling probe positions is shown in Figure 51. The 
concentrations of CO are approximately the same for the 90° and 
45° probe positions, except at distances between 0.08 cm (0.03 in.) 
and 0.23 cm (0.09 in.) from the surface where the CO concentration is 
lower for the 90° probe position. The CO concentration at the 0° 
probe position is the highest. The data for the 45° and 90° posi- 
tions are plotted in Figures 52(a) and (b). If the data are extrapo- 
lated to the surface, it can be seen that the surface concentration 
for the 45° probe position is greater than that at the 90° posit- 
ion. Therefore, the CO concentration at the surface increases with 
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position from the 90° to 0°. The increase in the CO concentration 
with decreasing angle can be attributed to convection of the combust- 
ion products upward and to the C0 2 reduction reaction with carbon, 
especially at the 0° probe position. 

The data from Figures 50 to 52 indicate that free convection 
played a major role in the combustion process. The decrease of 0 2 
with decreasing angle indicates that convected 0 2 was reacting at 
the surface and that the 0 2 could not be replenished by molecular 
diffusion. The increase in CO^ with decreasing angle also shows 
this same result, even though it was reacting at the surface to pro- 
duce CO. The data for the 0 2 indicates a change in control mechan- 
ism from chemical react ■ n rate control at the lower sampling posi- 
tions to mass transfer control at the top of the rod. This change of 
mechanism may also occur for the C0 2 at the 0° position, but data 
near the surface that would be needed to verify this, could not be 
obtained. 

Four tests were conducted in which a monochromatic optical pyrom- 
eter was used to determine the surface temperature. No gas samples 
were taken during these tests because the carbon rod holder assembly 
had to be removed from the mass spectrometer system. The same flow 
system was used and the conditions were kept the same as in the sam- 
pling tests. The time average surface temperature of the rod was found 
to be 1269* C, which is above the minimum surface temperature found in 
the literature for a gas phase reaction. There was only a few degree 
temperature difference between the bottom of the rod and the top of 


the- rod. 
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E. Comparison with Previous Experimental Data 

The concentration profiles obtained above for the horizontal rod 
can be compared to the concentration data of Wlcke and Wurzbacher [5] 
for the midpoint of a 6 cm vertical rod with a surface temperature of 
1230* C. In Figures 53 and 54 for the 90° and 45° probe posi- 
tions, respectively, It can be seen that the C0 2 concentrations are 
considerably higher and the 0 2 concentrations considerably lower for 
the vertical rod than for the horizontal rod. These differences are 
attributed to convective mass transfer. When compared to a diameter 
of 0.615 cm for the horizontal rod, the 6 cm length of the vertical 
rod allowed a greater distance for the build up of combustion pro- 
ducts, especially C0 2 . This would cause the increase in the COg 
concentrations and the decrease in the 0 2 concentrations. The CO 
concentrations for the vertical rod are higher than that for the hori- 
zontal rod, but the difference is considerably less than that for the 

C0 2 and O 2 concentrations. This indicates that CO is least 
affected by convective mass transfer. In Figure 55, the data taken at 

the 0° probe position for the horizontal rod are compared to the 
data for the vertical rod. Since convective effects for the horizon- 
tal rod were maximized at this 0° position, the data for the two 
rods should be more in agreement than for the 90° and 45° probe 
positions. This is seen to be true in Figure 55. For the vertical 
rod the C0 2 concentration rises more rapidly from the surface, 
reaching a maximum at approximately 0.25 cm from the surface, while 

the C0 2 concentration for the horizontal rod reaches a maximum at 
approximately 0.36 cm from the surface. For the vertical rod the 
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Figure 53. - Comparison of horizontal rod data with ver 
tical rod data * dtf* probe position. 





Figure 54. - Comparison of horizontal rod data with ver- 
tical rod data - 45° probe position. 




CO 



C0 2 concentrations were higher and the CO concentrations lower than 
fcr the horizontal rod. The 0 2 concentrations for the vertical rod 
are initially lower than that for the horizontal rod, but become high- 
er at approximately 0.36 cm from the surface. These differences in 
the data can be caused by the differences in the boundary layers for 
the two rod orientations. The boundary layer above the horizontal rod 
should be considerably thicker than the boundary layer for the verti- 
cal rod. This would cause the CO2 maximum to shift to a greater 
distance from the surface, the CO and C0 2 concentrations to decrease 
less rapidly, and the 0 2 concentration to increase less rapidly with 
distance for the horizontal rod. The large decrease in C0 2 , the 
decrease in CO, and the large increase in 0 2 for the horizontal rod 
at distances greater than 0.64 cm from the surface is attributed to 
ambient 0 2 mixing with the combustion gases. In general, for both 
rod orientations, the C0 2 concentrations increase with distance from 
the surface, reach a maximum and then decrease. Correspondingly, the 

0 2 concentrations increases with distance, and the CO concentrations 
decrease with distance from the surface. Considering the geometrical 
differences between the two experiments, the comparison of the data 
shows very good agreement as to the phenomena occurring during the 
combustion process. 

F. Conclusions for Normal Gravity Tests 

For carbon rods burning in a dry oxygen environment, there was a 
gas phase reaction in which CO was burned to C0 2 . From the 90° 
probe position to the 0° (top of the rod) probe position, the 0 ? 
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concentration decreased, the C (>2 concentration increased, but is 
expected to be lowest near the surface at the top of the rod, and the 
CO concentration, in general, increased. This means that convective 
effects were playing a considerable ole in the combustion process, 
decreasing the oxygen supply to tn. ■ A by convecting up the combus- 
tion products. This convective mass transfer process was considerably 
faster than the radial mass transfer process. 

From the data presented in Table 2, it can be concluded that for 
a 0.615 cm diameter rod with a surface temperature of 1269° C, the 
combustion process should be controlled by the mass transfer process. 
However, the normal gravity data obtained in this study indicates that 
the combustion of the carbon rod was controlled by the chemical pro- 
cess, except near the top of the rod where it was controlled by mass 
transfer. The change in control mechanism is caused, as mentioned 
above, by the combustion products being convected upward so the O 2 
could not easily reach the surface at the 0° probe position. It 
should be mentioned that previous investigators have used the carbon 
burning rate to determine the controlling mechanism, and did not use 
gas sampling techniques. 

Good agreement, especially at the 0° probe position, was ob- 
tained between the experimental data from this study and that of Wicke 
and Wurzbacher [5], showing that both studies have observed the same 
burning phenomena. 


V. MATHEMATICAL MODEL 


In order to compare the normal gravity data with a high temper- 
ature mathematical model that accounts for a gas phase reaction, the 
analysis of Caram and Amundson [6] for slabs and spheres was extended 
to cylindrical rods. The following assumptions were made: (1) a 

stagnant film of thickness 6 surrounds the carbon rod; (2) axial end 
effects are neglected; (3) the rod does not react internally and is 
impervious to mass transfer; (4) the species present are CO, C0 2 , 
and 0 ? ; (5) gas properties are constant; (6) the Schwab-Zeldovich 
approximation (that is, the diffusion coefficients are the same and 
equal to U/pc p )) holds; (7) the process is steady state; (8) ef- 
fective binary diffusion; (9) the flow is only in the radial direct- 
ion; and (10) there are two surface reactions, 2C + 0 2 » 2C0, and 
C + C0 2 » 2C0, and a gas phase reaction CO + 1/2 0 2 » C0 2 . The 
model is shown in Figure 56. The complete mathematical model deriva- 
tion is given in the Appendix, with only the essentials of the solu- 
tion being presented here. 

A. Governing Equations and Numerical Solution 


The equations governing this combustion process are as follows: 



Mixture Continuity Equation: 


prv - constant « m - M c r $ (2^ +^ co ) (1) 

2S 2S 


Species Continuity Equations: 


Oxygen: 



r dr 




2 M 0 2 **C0 


Carbon Dioxide: 


. dY C0, 
m w 2 

r Hr 


m d V‘0 

Carbon Monoxide: ~ —^r 



and the Thermal Energy Equation: 

!!!ffilb dl . ilb d_ / dr\ H C0^0 
r dr r dr y dry 2 


where p 
r 
v 


= density of the mixture; 

= radial distance from the center of the rod 
= radial velocity; 
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m 


M, 


Mr 


M 


M, 


CO r 

t 

CO 


2S 

*co 2S 

\ 

R co, 


'CO 


a, 


CO 


CO. 

f co 


CO. 

c 


CO 


■ carbon burning rate 

■ mole.ular weight of carbon; 

= molecular weight r)f oxygen; 

« molecular weight of carbon dioxide; 

* molecular weight of carbon monoxide; 

* surface molar reaction rate of oxygen; 

■ surface molar reaction rate of carbon dioxide; 

* gas phase mass reaction rate of oxygen; 

* gas phase mass reaction rate of carbon dioxide; 

= gas phase mass reaction rate of carbon monoxide; 

« gas phase molar reaction rate of carbon monoxide; 

* mass fraction of oxygen; 

= mass fraction of carbon dioxide; 

= mass fraction of carbon monoxide; 

» effective binary diffusivity of oxygen; 

= effective binary diffusivity of carbon dioxide; 

= effective binary diffusivity of carbon monoxide; 

= specific heat of mixture; 

= temperature at the edge of the stagnant film; 

= dimensionless temperature T /T b ; 

« thermal conductivity; and 

= heat of combustion of two moles of carbon 
monoxide. 


The surface boundary conditions at r * r^ assume that the flux 
of a species at the rod surface is equal to the surface reaction 
rate. The surface boundary conditions are 
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Oxygen: - 


pD o 2 dv o 2 
1 CTIT 


raY 0 

S + *0^"^2S 


( 6 ) 


p 0 co 2 dY co 2 

Carbon Oioxide: - 

co 2 


" Y co s 

S M CO, r S "^ C0 2S 


(7) 


p^Cq dY CQ 

Carbon Monoxide: — n j — 

M C0 dr 


mY CO 

+ JTir • 2 
S M CO r S 


(^2S + ^°2s) 


( 8 ) 


Thermal Energy: -XT fa £ - - H cofcO zs 


(9) 


where h 


0, 


TO, 


= heat produced by one mole of oxygen reacting with 
carbon to produce two moles of carbon monoxide; 
and 

* heat absorbed by one mole of carbon dioxide 
reacting with carbon to produce two moles of 
carbon monoxide. 


The boundary conditions at the edge of the stagnant film at 


r * r fa are: 


Oxygen: Y n = Y n 

U 2 °2b 


( 10 ) 


8 Y s 1 . V 

C0 2 C0 2b °2b 


Carbon Dioxide: Y 


(ID 
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Carbon Monoxide: • 0 

(12) 

Thermal Energy: t ■ 1 

(13) 


where Y n and Y rn are the ambient concentrations of oxygen 
u 2b LU 2b 

and carbon dioxide, respectively. Equations (1) through (13) give a 
complete mathematical formulation of the problem. 

Equations (2) and (3), (4) and (3), and (5) and (3) were combined 
to eliminate The equations were solved yielding 



■ 2bt * t », * K* • ” * c * sJfey • 


20t$ - e 


(16) 


where 


_ m 

Y cl! 


M co ? c P T b 

0 a n- ■ ■ 

H co 
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If one of the variables V Q . Y™ . V 


. -eg . - co . or t is known, 


2 2 

the other variables can be determined using equations (14), (IS), and 
(16). In this case, the surface temperature was considered known, and 


the surface mass fractions Y 


o 2S - V C0 2S - 


and Y C g were de- 


termined. The gas phase reaction between carbon monoxide and oxygen 
was assumed to be first order with respect to both carbon monoxide and 
oxygen. 


% " k g C C0 C 0 2 e *P(- E / RT ) 


(17) 


where kg 



£ 


R 


« pre-exponential factor for the gas phase reaction, 
CO ♦ 1/2 0 2 » C0 2 ; 

« molar concentration of carbon monoxide; 

» molar concentration of oxygen; 

» activation energy for the gas phase reaction, 

CO + 1/2 » C0 2 ; and 

« universal gas constant. 


This equation was substituted into equation (4) yielding 


m 



pTT' 



- M C0 k g C C0 C 0 2 exp(-E/RT) 


(18) 


Equation (Iff) was solved using a Runge-Kutta technique. The solu- 
tion of equation (18) along with equations (14), (15), and (16) yield- 
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ed carbon monoxide, carbon dioxide, and oxygen concentration profiles, 
and the temperature profile. To compare the results from the model 
with the normal gravity data, the mass fractions were converted to 
mole fractions using the equation 

X 1 ■ (W ♦ (Vj/Mj) ♦ (yn k ) I* 9 ) 


where i represents oxygen, carbon dioxide, or carbon monoxide, and j 
and k represents the other two species. 

The physical properties used in the equations are presented in 
Table 6. The binary diffusion coefficients used were obtained by av- 
eraging the coefficients obtained from the Stef an-toaxwel 1 equations 
for an ideal gas. The activation energy for the gas phase reaction of 
carbon monoxide with oxygen was taken as 43800 cal /g mole [40]. The 
rate constants for the surface reactions, that is, for carbon with 
oxygen and carbon with carbon dioxide, were taken from Reference [41] 
and are listed in Table 7. Only the reaction rate constants that 
allowed physically reasonable solutions, that is, surface mass frac- 
tions equal to or greater than zero and less than or equal to one us- 
ing equations (14), (15), and (16) are listed. 

In the solution of the equations using the Runge-Kutta technique, 
the pre-exponential factor, k ? , for the gas phase reaction was var- 
ied so that the boundary conditions at the outer edge of the stagnant 

film were met. These were initially taken to be Y n . i 

u 2b 

Y rn * 0, and r • 1. The thickness of the stagnant film was var- 
LU 2b 

ied from 0.302 cm (0.119 in.) to 0.759 cm (0.299 in.) to account for 


Table 6. - Physical properties 


f 

c p 

D 

T b 

«co 

s 

»co. 

T a 


4, 982x10“* g/cra^ 
0,273 cal/g °K 
1,23 cm 2 /sec 
294,44 °K 

!• 348x10^ cal/g mole 
2,677x10* cal/g mole 
4,063x10* cal/g mole 


« 1542 
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Table 7. - Surface reaction rate 


conetante 


k ft • 1.1 x 10^ exp(-25500/RT) cm/sec 
°2 

k n * 1,0 x 10^ exp(-20300/RT) cm/e'c 
°2 

k * 6.9 x 10 6 exp(-44250/RT) cm/eec 

co 2 

k C Q « 1.6 x 10^ exp(-52000/RT) cm/sec 

2 

k C o * 3*0 x 10 6 exp(-40400/RT) cm/sec 
k co = 3.16 x 10 10 exp ( -74000/RT ) cm/sec 
38 7.9 x 10^ exp(-51200/RT) cm/'sec 
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the Increase In the free convective boundary layer around the circum- 
ference of the horizontal carbon rod. The stagnant film thickness and 
surface reaction rate constants used in the solution of the model 
equations, along with the resulting surface mole and mass fractions of 
oxygen, carbon dioxide, and carbon monoxide, the derivative of the 
mass fraction of CO at the surface, and the pre-exponential factor for 
the gas phase reaction are presented in Table 8. The computed mole 
fractions were compared to the normal gravity experimental data ob- 
tained at the 80° and 45° probe positions. The 0° p 'be posi- 
tion data was not used for comparison, since it was not possible to 
obtain concentrations close to the surface at that probe position. 

B. Comparison of Model with Normal Gravity Experimental Data 

The cases where the model compared best with the normal gravity 
experimental dcta are presented in Figures 57 to 60. In Figures 57 
and 58, the experimental data taken at the 90° and 45° probe posi- 
tions, respectively, are compared to the mole fractions calculated 
using a stagnant layer thickness of 0.404 cm (0.159 in.). In Figure 
57, it can be seen that the model does predict the mole fraction be- 
havior of CO fairly well for the 90° probe position, especially at 
distances greater than 0.071 cm (0.028 in.) from the surface. The 
model over predicts the mole fractions of 0 ^ and under predicts the 
mole fractions of COg. For the 45° probe position, shown in Fig- 
ure 58, the calculated mole fractions of CO are higher than the exper- 
imentally obtained mole fractions at distances less than 0.076 cm 
(0.03 in.) from the surface, and are lower than the experimental mole 


T*bU 9# - flodol parwi»t«r« «,i4 
rsaultc 


o.io? 

(o.m) 

7-1 

0. 102 

(0.119) 

7-1 

0. 10? 

(0.U9) 

7-1 

0. )02 

(0.119) 

7-? 

0. 10? 

(0.119) 

7-2 

0.151 

(0.1)9) 

7-1 

0. 151 

(0.1)9) 

7-1 

0.153 

(0.1)9) 

7-1 

0. 151 

(0.1)9) 

7-2 

0. }51 

(0.1)9) 

7-2 

0.404 

(0.159) 

7-1 

0.404 

(0.159) 

7-1 

0.404 

(0.159) 

7-1 

0.404 

(0.159) 

7-2 

0.404 

(0.159) 

7-2 

0.455 

(0.179) 

7-1 

0.455 

(0.179) 

7-1 

0.455 

(0.179) 

7-1 

0.455 

(0.179) 

7-2 

0.455 

(0.179) 

7-2 

0.501 

(0.198) 1 

1 7-1 

0.503 

(0.198) 

7-1 

0.50) 

(0.198) 

; 7-1 

0.50) 

(0.190) 

I 7-2 

0.503 

(0.198) 

7-2 


m 


1.48x10 

4.37x10 

5.74x10 


•urfaoa raaotlon rat* constant llatad In Tabla 7 















122 


Tabl* 6*. Concluded 


i 

—at 

Jo, 

*co 2 

St 

Hs 

*co s 

*0, 

2< 

• x *». 

. 

■ 

0 

0 

X 

*W.r 

1/on (1/lnJ 

y 

nol« »*c 

0.556 

(0.219) 

7-1 ’ 

7-3 * 

0.769 

0.150 

0,081 

0.79? 

0.112 

0.095 

-o.aj9 

(-2.131) 

V.076 

2. 4 3x10* 3 

0.556 

(0.219) 

7-1 

7-4 

0.537 

0.039 

0.424 

0.51? 

0.027 

0.461 

-2.646 

(-6.727) 

0.2069 

1.96X10 11 

0.556 

(0.219) 

7-1 

7-5 

0.743 

0.138 

0,199 

0.694 

0.094 

0.212 

-1.045 

(-2.654) 

0.090 

1.37x10* 3 

0.556 

(0.219) 

7-2 

7-6 

0.618 

0.078 

0.303 

0.605 

0.056 

0.339 

-1 .088 

(-2.764) 

0.159 

3.97x1^ 

0.556 

(0.219) 

7-2 

7-7 

0.6002 

0.0694 

0.3305 

0.5837 

0.0491 

0.3672 

-1.258 

(-3.196) 

0.1695 

4.88.10* 2 

0.607 

(0.239) 

7-1 

7-3 

0.763 

0.147 

0.090 

0.784 

o.no 

0.106 


(-2.091) 

0.075 

1.91xl0* 2 

0.607 

(0.239) 

7-1 

7-4 

0.5328 

0.037) 

0.4)00 

0.5069 

0.0258 

0.4673 

-2. * ?? 

(-6.419) 

0.1984 

1.57x10** 

0.607 

(0.239) 

7-1 

7-5 

0.736 

0.135 

0.129 

0.750 

0.100 

0.150 

-1.022 

(-2.597) 

0.089 

1.03K10 13 

0.607 

(0.239) 

7-2 

7-6 

0.609 

0.073 

0.318 

0.594 

0.052 

0.354 

—1,061 

(-2.696) 

0.156 

3.45xlO* 2 

0.607 

(0.239) 

7-2 

7-7 

0.591 

0,065 

0.345 

0.572 

0.046 

0.332 

-i.r? 

(-3.092) 

0.166 

4.12*10 1? 

0.656 

(0.259) 

7-1 

7-3 

0.757 

0.145 

0,098 

0.777 

0.108 

0.115 

-0.812 

(-2.062) 

0.074 

1.56x10* 3 

0.656 

(0.259) 

7-1 

7-4 

0.5294 

0.0357 

0.4350 

0.5031 

0.0247 

0.4723 

-2.422 

(-6.151) 

0.1909 

1.30x10** 

0.656 

(0.259) 

7-1 

7-5 

0.731 

0.132 

0.137 

0.74 3 

0.096 

0.159 

-1.002 

(-2.545) 

0.007 

9.25**.2 1 3 

0.656 

(0.259) 

7-2 

7-6 

0.596 

0.069 

0.333 

0.581 

0.049 

0,370 

-1.031 

(-2.619) 

0.154 

2.75xl0* 2 

0.658 

(0.259) 

7-2 

7-7 

0.582 

0.061 

0.357 

0.563 

0.043 

0.394 

-1.180 

(-2.996) 

0.163 

3.14xlO* 2 

0.709 

(0.279) 

7-1 

7-3 

0.753 

0.142 

0.105 

0.771 

0.106 

0.123 

-0.798 

(-2.026) 

0.073 

1.29xl0 13 

0.709 

(0,279) 

7-1 

7-4 

0,5263 

0.0343 

0.4395 

0.4996 

0.0237 

0,4767 

-2.330 

(-5.917) 

0.1843 

l.03*10 11 

0.709 

(0.279) 

7-1 

7-5 

0.726 

0.129 

0.145 

0,737 

0.095 

0,168 

-0.901 

(-2.491) 

0.086 

7.37*10 12 

0.709 

(0.279) 

7-2 

7-6 

0.591 

0.065 

0.344 

0.573 

1 0.046 

0.351 

-1.012 

(-2.570) 

0.151 

2.51xl0 12 

0.709 

(0.279) 

7-2 

7-7 

0.575 

0.057 

0,368 

0.555 

0.040 

0.405 

-1.147 

(-2.914) 

0.159 

3.29xlO* 2 

0.759 

(0.299) 

7-1 

7-3 

0.746 

0.140 

0.11? 

0.765 

1 0,104 

0.131 

-0.757 

(-1.998) 

0.073 

9.00xl0* 2 

0.759 

(0.299) 

7-1 

7-4 

0.524 

0.033 

0.443 

0.497 

0.023 

0.480 

-2.244 

(-5.701) 

0.178 

7.95x10*° 

0.759 

(0.299) 

7-1 

7-5 

0.720 

0.127 

0.153 

0,729 

0.094 

0.177 

-0.965 

(-2.45?) 

0,085 

6.29xl0 12 
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(0.299) 

7-2 

7-6 

0.503 

0.061 

0.356 

0.564 

0.043 

0.393 

-0.989 

(-2.512) 

0.149 

2. 13x10* 2 
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(0,299) 

7-2 

7-7 

0.5674 

0.0537 

0.379 

0.5459 

0.0376 

0.4166 

-1.115 

(-2.833) 
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ko 2 • L U10 3 «xp(-2550CKRT) cm/sec 
k CO, m 1.6*10 9 e«p|-5200QIRTI cm/sec 
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Figure 57. - Comparison of model with experimental data - 0. 405 cm 
<0. 159 in. ) stagnant film thickness - probe position 90°. 



TEMPERATURE, K 
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‘ Comparison of model with experimental data - a 404 
iO. 159 in. ) stagnant film thickness * probe position 45° 
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fractions at greater distances. Again, the model over predicts the 
mole fractions of and under predicts the mole fractions of CO 2 . 

In Figures 59 and 60, the experimental data taken at the 90° and 
45° probe positions are compared to the mole fractions calculated 
using a larger stagnant layer thickness of 0.503 cm (0.198 In.). 

Figure 59 shows that at the the 90° probe position, the model pre- 
dicts CO mole fractions slightly higher than those obtained experimen- 
tally. At the 45° probe position, shown In Figure 60, the model 
predicts higher CO mole fractions than obtained experimentally at dis- 
tances less than 0.107 cm (0.042 in.) from the surface, between dis- 
tances of 0.107 cm (0.042 in.) and 0.203 cm (0.080 in.), there is good 
agreement in the mole fractions, and at distances greater than 0.203 
cm (0.080 in.) the model predicts slightly higher CO mole fractions. 

At both probe positions the model over-predicts the 0£ mole frac- 
tions and under-predicts the C0 2 mole fractions. 

In an attempt to obtain better agreement between the experimental 
data and the model, the mole fractions of 0 2 and C0 2 at the outer 
edge of the stagnant film were changed to 0.925 and 0.075, respective- 
ly. The values were reasonable because of the design of the burning 
apparatus, which could have caused some recirculation of the combust- 
ion gases and thereby increased the mole fraction of CO? in the am- 
bient atmosphere. In the solution of the equations only the surface 

3 

reaction rate constants k n = 1.1 x 10 cm/ sec exp [-25500/RT] 

u 2 

and k^ 0 * 1.6 x 10 9 cm/ sec exp [-52000/RT] were used because 
these gave the closest agreement for the case of an oxygen mole frac- 
tion of 1.0 at the outer edge of the stagnant film. The thickness of 
the stagnant film was again varied from 0.302 cm (0.119 in.) to 0.709 
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cm (0.279 in.). The stagnant film thickness and the surface reaction 
rate constants used in the solution of th'.< model equations, along with 
the resulting surface mole and mass fractions of oxygen, carbon di- 
oxide, and carbon monoxide, the derivative of the mass fraction of CO 
at the surface, and the pre-exponential factor for the gas phase re- 
action are presented in Table 9. The cases where the model and the 
experimental data compare the best are presented in Figures 61 to 64. 

In Figures 61 and 62, the experimental data taken at the 90° and 
45° probe positions are compared to the mole fractions calculated 
using a stagnant layer thickness of 0.353 cm (0.139 in.). Figure 61 
shows that at the 90° probe position, the calculated mole fractions 
of CO and CO^ and the experimental data are in good agreement at 
distances greater than 0.117 cm (0.046 in.) from the surface. The 
calculated and experimentally obtained 0^ mole fractions are in fair 
agreement over the whole distance experimental data was taken. At the 
45° probe position, shown in Figure 62, the model predicts 0^ mole 
fractions that are higher and CO and CO^ mole fractions that are 
lower than the experimental ly obained mole fractions. 

In Figures 63 and 64, the data obtained at the 90° and 45° 
probe positions are compared to the mole fractions calculated using a 
stagnant layer thickness of 0.429 cm (0.169 in.). For the 90° probe 
position, shown in Figure 63, the model predicts 0 ? mole fractions 
that are lower and CO and CO^ mole fractions that are higher than 
the experimentally obtained mole fractions. Figure 64 shows that at 
the 45° probe position, the calculated 0£ mole fractions are in 
good agreement with the experimentally obtained data at distances 
greater than approximately 0.051 cm (0.02 in.) from the surface. The 
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Figure 61. - Comparison of model with experimental data - 0. 353 cm 
(0. 139 in. ) stagnant film thickness - probe position 9flP. 
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*0 2 ■ 1. IxKpexpl-SSOQfRTI cm/sec 

k C0 2 • L 6xl0 9 exp(-MOO0RTl cm/sec 

l • 4X00 cal/g molt 
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Flgurt 64. - Comparison ol model with experimental data - ft 429 cm 
(0. 169 in. ) stagnant film thickness - probe position 45°. 
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calculated mole fractions of CO and CO2 are in good agreement with 
the experimental data for distances greater than approximately 0.140 
cm (0.055 in.) from the surface. 

It is consistent with the boundary layer concept that the boundary 
layer thickness should increase from the 90° to the 45° probe 
positions. This fact is reflected in Figures 61 and 64, which show 
that the 90° data is best fit with a film thickness of 0.353 cm, 
while the 45° data is best fit with a film thickness of 0.429 cm. 

To summarize the model comparison with the experimental data, the 
model with an 0^ mole fraction equal to 1 at the outer edge of the 
stagnant layer predicted the mole fractions of CO over the distance 
from the surface better than the model with mole fractions of O2 
equal to 0.925 and CO2 equal to 0.075 at the outer edge of the stag- 
nant film. However, it did not accurately predict the mole fractions 
of CO2 nd O2. The model with mole fractions of O2 equal to 
0.925 and CO2 equal to 0.075 predicted the mole fractions of O2 
reasonably well vp t 'e distance from the surface where data was 
taken. Also, it predicted the CO2 and CO mole fractions at dis- 
tances greater than approximately 0.129 cm (0.051 in.) from the sur- 
face. Using either outer edge boundary condition, the model was un- 
able to predict the experimentally obtained C0 2 mole fractions at 
distances less than approximately 0.1 cm (0.039 in.) from the sur- 
face. This is most probably due to free convective effects, which 
caused the combustion products to rise. The experimental data taken 
at the 90° and 45° probe positions show that the largest differ- 
ence in both the O2 and CO2 mole fractions occurs near the surface 
and decreases with distance from the surface. This means free convec- 
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tion mass and heat transfer effects would be greatest near the surface 
and would decrease with distance from the surface. The stagnant layer 
thickness was varied in the model in order to somewhat account for the 
convective effects, and did improve the model agreement with the 
data. The model thus qualitatively describes the carbon rod combus- 
tion process, and should either be compared with zero-gravity combus- 
tion data or modified to more fully include convective effects. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


Experiments were conducted in zero and normal gravity in which 
horizontally mounted spectroscopic carbon rods were burned in dry oxy- 
gen environments. 

In the zero gravity experiments 0.615 cm (0.242 in.), 0.457 cm 
(0.180 in.), and 0.305 cm (0.120 in.) diameter carbon rods were burn- 
ed in oxygen at pressures of 3.45 x 10 4 N/M 2 (5 psia), 6.89 x 
10 4 N/M 2 (10 psia), 1.03 x 10 5 N/M 2 (15 psia), and 1.38 

5 p 

x lCr N/M £ (20 psia). The oxygen had a maximum water vapor 
content of 0.5 ppm. The photographs of the zero-gravity combustion 
process showed that in all the cases where sustained burning occurred, 
a blue flame extending from the surface surrounded the rod. This 
means that a gas phase reaction was taking place in which carbon mon- 
oxide was being oxidized to carbon dioxide. Mathematical models des- 
cribing the carbon combustion process which account for a gas phase 
reaction are those of Hugo, Wicke, and Wurzbacher [4] and Caram and 
Amundson [6]. It was not possible to test these models since concen- 
tration profiles and surface temperature measurements could not be 
obtained in the zero gravity tests. From the zero gravity data, it 
was found that the ratio of flame diameter to rod diameter, d*, 
increased as the pressure decreased, and decreased as the rod diameter 
increased. The data was used to obtain the correlation equation: 
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d* = o(l - e"*l t_ *2) 

where a « 1.943 - 6.599(d - 0.18) - 0.03517 (P - 12.5) ♦ 
46.69(d - 0.18) 2 
- 0.687 

0 2 - 1.015 ♦ 4. 335(d - 0.18) * 0.0278(P - 12.5) - 
119. 4(d - 0.18) 2 


and d (rod diameter) is in inches, P (pressure) is in psia, and t 
(time) is in seconds. The zero gravity data obtained in this study 
can be used to properly design a long duration Spacelab experiment on 
carbon combustion, in which complete information on concentrations and 
temperatures can be obtained. 

In the normal gravity experiments, 0.615 cm (0.242 in.) diameter 
carbon rods were burned in oxygen at atmospheric pressure. The oxygen 
had a moisture content of less than 6 ppm. A concentration sampling 
probe was positioned near the circumference of the rod at either the 
top or at angles of 45° or 90° from the top, and yielded concen- 
tration profiles of C0 2 , CO, and 0 2 as a function of distance from 
the surface. At all the sampling probe positions, the 0 2 concentra- 
tion increased with distance from the surface, while the CO concentra- 
tion decreased with distance. At the 90° position the C0 2 
concentration first increased slightly with distance from the surface 
and then decreased, whereas at the 45° probe position, it decreased 
with distance from the surface. At the 0° probe postion, the C0 2 
concentration increased to a maximum away from the surface, and then 
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decreased. Moving from the 90° probe position to the top probe 
position, the 0^ concentration in the gas surrounding the rod de- 
creased, the C0 2 concentration increased, but is expected to be low- 
est near the surface at the top of the rod, and the CO concentration 
in general increased. This meant that convective effects were playing 
a considerable role in the combustion process, decreasing the oxygen 
supply to the rod by convecting up the combustion products. The data 
also indicates that the surface reactions at 90° and 45° from the 
top of the rod were controlled by the chemical process, while the sur- 
face reaction at the top of the rod may have been controlled by mass 
transfer. The data shows that there was a gas phase reaction in which 
CO was burned to C0 2 , as was the case in the zero-gravity experi- 
ments. The time averaged surface temperature of the rod was found to 
be 1269° C, which is above the minimum temperature found in the lit- 
erature for a gas phase reaction to occur. 

The carbon dioxide mole fraction increasing with distance from the 
surface and then decreasing at the 90° and 0° probe positions 
means that C0 2 was being reduced to CO at the surface. This re- 
duction reaction was probably occurring at the 45° probe position 
also, but was masked by free convective effects which caused the com- 
bustion products, especially C0 2 , to rise. The C0 2 reduction re- 
action occurring in these experiments with a carbon surface temper- 
ature of 1269° C is in agreement with the finding of other 
investigators, who determined the minimum temperaure at which this 
reaction would occur to be 1100°C. 

Good qualitative agreement between the carbon combustion data ob- 
tained in the present investigation at the 90° and 45° probe posi- 
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tions and that obtained experimentally by Wicke and Wurzbacher [5] 
indicates that the same physical phenomena was being observed in both 
studies. This conclusion is justified by quantitative agreement of 
the data at the 0° probe position, where convective effects were of 
similar importance in both studies. 

The normal gravity experimental data were also compared to the 
stagnant film mathematical model proposed by Caram and Amundson [6]. 
The model, with ambient mole fractions of 0^ equal to 0.925 and 
C0 2 equal to 0.075, predicted the mole fractions of 0 2 reasonably 
well over the complete distance from the rod surface, and the mole 
fractions of C0 2 and CO at distances greater than approximately 
0.129 cm (0.051 in.) from the surface. The stagnant film thickness 
was varied in the model in order to somewhat account for the convect- 
ive effects, and did improve the model agreement with the data. The 
model thus qualitatively describes the carbon rod combustion process, 
and should either be compared with more complete zero-gravity data or 
modified to more fully include convective effects. 

It is therefure concluded, in agreement with other investigators, 
that carbon combustion is complex. Under the experimental conditions 
used in the present study, the combustion process was characterized by 
two surface reactions, 2C + 0 2 » 2C0 and C + C0 2 » 2C0, and a gas 
phase reaction, CO + 1/2 0 2 » C0 2 . The process controlling step 
changed from chemical process control to mass transfer control as 
probe measurements went from 90° to 0° in the normal gravity ex- 
periments. In both zero and normal gravity, the carbon combustion is 
best described mathematically by the model of Caram and Amundson [6], 
although modification for convective effects is needed. 
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To learn more about the role of convection as it affects the mech- 
anisms of carbon combustion, a long-term zero gravity experiment is 
needed. This would eliminate the strong normal gravity convective 
effects and allow for long term gas sampling and temperature measure- 
ments. The obtained surface temperature, concentration profiles, and 
temperature profiles would completely describe the zero-gravity com- 
bustion process and confirm the zero-gravity combustion model. By 
comparing normal gravity and zero-gravity data, the extent that free 
convection affects the Og, CC^, and CO concentrations can be de- 
termined. This information can be used to obtain a correction factor 
for the stagnant film model so that it could predict the concentration 
profiles when convection is present. Also, the zero-gravity data can 
be used to obtain information concerning the relationship between mass 
transfer and chemical process control. Additional supplementary ex- 
perimental work is needed to obtain accurate reaction rate constants 
for the carbon-oxygen and carbon-carbon dioxide surface reactions and 
for the oxygen-carbon monoxide gas phase reaction. 


APPENDIX 


MATHEMATICAL MODEL DERIVATION 


The model presented by Caram and Amundson [6] was extended to a 
cylindrical rod. The following assumptions were used: (1) a stagnant 

film of thickness 6 surrounds the carbon rod; (2) axial end effects 
are neglected; (3) the rod does not react internally and is impervious 
to mass transfer; (4) the species present are CO, CO^, and 0 2 ; (5) 
gas properties are constant; (6) the Schwab-Zeldovich approximation 
(that is, the diffusion coefficients are the same and equal to 
(A/pCp)) holds; (7) the process is at steady state; (8) effective 
binary diffusion; (9) the flow is only in the radial direction; and 
(10) there are two surface reactions, 2C + 0 2 * 2C0, and C +C0 2 » 

2C0, and a gas phase reaction, CO + 1/2 0^ » C0 2 . 

A. Governing Equations 


The equations governing this combustion process are the Species 
Continuity Equations: 


3p 


It * r tr (rn 1r ) * R 1 1 ■ 0 2 ,C0 2 ,C0 


(Al) 
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where 



(A2) 


and the Thermal Energy Equation: 



(A3) 


where p^ 
t 
r 



T 


x 


= mass concentration of species i; 

= time; 

= radial distance from the center of the rod; 

= mass flux of species i in the radial direction; 
* gas phase mass reaction rate of species i; 

= effective binary diffusivity; 

= radial velocity; 

= density of mixture; 

= specific heat of mixture; 

= temperature; 

= thermal conductivity of the mixture; 


q' 


'CO 


* 


CO 


H co * rn 

= g = the heat generated by the combustion of 

carbon monoxide; 

= heat of combustion of two moles of carbon 
monoxide; and 

a gas phase molar reaction rate of carbon monoxide. 


By using the assumptions of steady state, constant physical prop- 


erties, and binary diffusion, and combining equations (Al) and (A2), 
the Species Continuity and Thermal Energy Equations become 





(Ml 


°V 


dT 

dr 



♦ q' 


(A5) 


where * mass fraction of species i; and 

. effective binary diffusivity. 

At steady state conservation of mass for the mixture yields 
(d/dr)(rpv) = 0. Therefore prv = constant = m is the carbon burning 
rate. 

Using equation (A4) and the fact that prv * m is a constant, the 
individual Species Continuity Equations can be written as 


r dr 


pO r 


* 2 d 
T~dF 



(A6) 


m 

r 



dr 




(A7 


.-it 


aitox»a» 


unit- 
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m ^CO _ OUcO _1 / ^[co\ . . 
r dr r dr \ dr ) 


The Thermal Energy Equation, (A5), can be written as 


7* 37 ' 7 If ( r 3?) 


H C0^C0 


Because mass must be conserved, the gas phase mass reaction rates in 
equations (A6) to (A8) must add up to zero, that is 


_r o 2 ’ R C0 + r co 2 " 0 


(A10) 


From the gas phase reaction formula, CO ♦ 1/2 0 2 > C0 2 , it can be 
seen that for every mole of C0 2 formed, one mole of CO and 1/2 mole 
of 0 2 are consumed. Therefore, in molar units, 

.1 » m ca (A1 1 ) 


(A12) 


The mass reaction rates, R^, and the molar reaction rates, 
are related by the equation 


R, ■ 
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where « molecular weight of species i. 
Equations (All) and (A12) can thus be written as 

r co 2 “ m co 2 ^co 
^ M 0 2 *co 


(A13) 


(A14) 


and, in the same manner, R CQ can be written as 

R C0* M C0^C0 (A15) 

Substituting equations (A14) into (A6), (A13) into (A7), and (A15) 
into (A8), and substituting r * T/T^, where Tj, * temperature at 
the edge of the stagnant film, into (A9), the following equations are 
obtained. 


m 

r 


dr 



' 2 M 0 2 ^C0 


(A16) 


dY 


m 

r 


CO. 

dr 



dY 


CO, 


m co 2 ^co 


r 


dr 


(A17) 
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m dY C0 p0 C0 d / dY C<>\ M a, 

7 i ? — r d? \ r it; ■ - M co^co 


A 


mc p T b dT xT b d / dt\ . H C0 ^ 

r dr r dr y dr/ ? 


CO 


B. Determination of the Burning Rate. 


(A18) 


(A19) 


To determine the burning rate of carbon, the surface reactions 

must be considered. From Figure Al, it can be seen that for every 

molecule of 0 2 that reacts at the surface, two molecules of CO are 

formed. Therefore, two atoms of carbon are removed from the surface. 

The molar rate of formation of CO at the surface is equal to twice the 

molar consumption rate of 0 2 at the surface. Therefore, the mos.. of 

carbon removed from the surface by 0, is M_2# n . From Fig- 

c c u 2s 

ure A2, it can be seen that for every molecule of C0£ that reacts at 
the surface, two molecules of CO are formed. But, because one carbon 
atom is carried to the surface by the C0 2 molecule, only one carbon 
atom is removed from the surface in forming two CG molecules. There- 
fore, the mass of carbon removed form the surface by CO2 is 


M # • The total mass flow from the surface Is, therefore, 

c co 2s 


* m c( 2 *o 2S * *C0 2S ) 
or 

m = M~r c |2^f n ♦ S rA j 
C % \ °2S C0 2S/ 


(A20) 


(A21) 



where = molecular weight of carbon; 

# n - molar surface reaction rate of oxygen; 

U 2S 

^COor - molar surface reaction rate of carbon dioxide; and 
2 S 

r<. = radius of carbon rod. 

C. Determination of Boundary Conditions. 


At the surface of the rod, it is assumed 

cules that hit the surface react. Therefore, 

flux N C() , 
u 2S 


that all the COg mole- 
in terms of the molar 




(A22) 


Because the molar flux and the mass flux are related by 


Mj 


9 


equation (A 22) can be written as 

n co 2S = ^CO^CO^ 


(A23) 


Using equation (A2), equation (A23) can be written as 


dY 


CO, 


-pD, 


C0 2 dr 


* pV C0 2s v = - m co/co 2S 


Using prv = m this equation becomes 


pD co 2 dY co 2 


M, 


CO, 


dr 


mY c° 2 s 

r S M C0. = ‘ ^ C0 2S 


(A24) 
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All the Op molecules that reach the surface react, so 


■% " *0 
°2S °2S 


(A25) 


By following the same procedure It can be shown that 


M n dr 
°2 


mY 0 

+ _I2S . 

s r S M 0 2 °2S 


(A26) 


Because every molecule of 0^ that reaches the surface forms two 
molecules of CO, and every molecule of COp that reaches the surface 
forms two CO molecules, the molar reaction rates are related by the 
equation 




CO, 


( s 


21*0 + ^CO 

U 2S C0 2S 


) 


Therefore, at the surface of the carbon, the flux of CO is 


N, 


CO, 


+ 2 K + '“») 


or 


n C0 s " 2M C0 


(%s * ‘Ss) 


(A27) 


Using equation (A2) and prv = m, equation (A27) can be written as 


p0 C0 dY C0 
M C0 dr 


mY 


CO, 


r S M C0 


* * <Hs) 


(A28) 


!/ m * « 


i 
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At the surface the heat flux is equal to the heat produced by the 


reaction, 2C ♦ O 2 » 2C0 (Hp^), and the heat absorbed by the re- 


action, C + CO 2 ■> 2C0 (H co ^), or 


x 41 

’ X df 


8 H o *0 " H co *co 

u 2 U 2S tu 2 lu 2S 


where H f . = heat produced by one mole of oxygen reacting with 

'2 

carbon to produce two moles of carbon monoxide; and 


H 


CQ * heat absorbed by one mole of carbon dioxide reacting 


with carbon to produce two moles of carbon monoxide. 
Using r = T/T^ this equation becomes 


■XT — 
A, b dr 


8 H o *0 
S U 2 U 2S 


H C0 2 *C0 2s 


(A29) 


At the outer edge of the stagnant film, the following conditions 
are assumed. At r = r. 


b> 


'CO 


Y o 

(A30) 

°2b 


Y C0 

tu 2b 

(A31) 

= 0 

(A32) 
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T - 1 


(A33) 


C. Solution of Covering Equations. 


The equations to be solved are (Alb), (A17), (A18), and (A19), and 
can be written as: 



(A16) 


(A17) 


m 

rM, 


dY 


CO 


CO 

dr 


*2 

rM, 


CO 

l C0 


JL L 

dr V dr / ^CO 


2mc o T b dT 
rH C0 dr 


2xT fc 


rH 


CO 


A L dl\ 

dr l r dr J 


r C0 


(A18) 


(A19) 


The boundary conditions at the surface are equations (A24), (A26), 
(A28), and (A29), at r = r $ ; 


pC C0 2 



mY 


CO 


2S 


r S M C0, 



(A24) 
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S 




(A26) 


p0 C0 dY C0 


M 


CO 


dr 


mY CO, 

+ nr* 
r s M co 



(A28) 


xT b dF 


* Hq “ ^co ^co 
S °2 °2S C0 2 C0 2S 


(A29) 


The boundary conditions at the outer edge of the stagnant layer are 


at r 


\ ' \ 

(A30) 

v co 2 ’ v co 2b 

(A31) 

-< 

o 

o 

It 

o 

(A32) 

T * 1 

{ A33) 


The equations for the boundary conditions were also combined. Equa- 


tion (A26) was multiplied by two and combined with equation (A24), 
equations (A28) and (A24) were combined, and equation (A29) was multi 
plied by 2 / h C q and combined with equation (A24) yielding 


"5 

iMwitr 
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2D« p dY« 2mY#» pD^a dY-» roY«« / 

°2 °2 . __2S C0 2 C0 2 . C0 2S _ /,* 

m o 2 s r s , \) 2 m co 2 ^ s r s M co 2 \ 1 


1 + Qt 

o 2S *co 2 < 


(A37) 


^co^!cg , co s 

M co dr r S M C0 

s 


pD C0 2 dV C0 2 roV C0 2s 
m co 2 dr s + r S M C0 2 


2^0 + ^rn 

°2S C0 2S 


2xT b d, C0c0 2 dYc0 2 "Hs 

H co ^ s * "co 2 dr s ‘ 


(A38) 


H cq \V°2S ' H C0 2 «C0 2S ; + *C0 2 


(A39) 


The boundary conditions at the outer edge of the stagnant film were 
likewise combined. The equations containing mass fractions were di- 
vided through by the respective molecular weights, and the equation 
containing r was multiplied through by - 2c p T b /H C Q. Also, the 

equation containing Y , that is, equation (A30), was multiplied 

u 2 

by two. Performing the above steps, then combining equations (A30) 
and (A31), (A31) and (A32), and (A31) and (A33) yields 


(A40) 


Vg °2 , Y C0 Sb 

M C0 2 *To H COj 


!2t . ifelfe t . !2tt . fffilb 

M C0 2 h co 1 m co 2 h co 


Equation (A34) was integrated, yielding 


2mY 0 2pD 0 dY 0 
°2 °2 _ °2 


°o 2 dY o 2 mV co 2 po co 2 dv co 2 

U r ~di “ + M co ' ” r dr * A 
°2 C0 2 M C0 2 


This equation, evaluated at the surface, gives 


2 ° V dV °2 + "Hs 


po co 2 dv co 2 

C0 2 S 


and, compared to the surface boundary condition, equation (A37) 
written as 


2 D o 2 p dY o 2 mY co 2i 

"V rs_ ^" s + 1 ^ 


r S dr 




2 %S * * C W S 


the value of the constant A is determined as 
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Using the Schwab-Zeldovich approximation Oq^ * * 0, 

equation (A43) can be written as 


2 l . _j- v t ^V\ t _!_ / v + dYc0 ? \ . . A 

M Q 2 l Y °2 dr J + M C0 2 l p®** Yc0 2 dr ] pTF 


Using the Integrating factor e* 
was Integrated, yielding 


/ Im/pDr)dr 


■ r' m ^ p ® , this equation 


t (v*’’) • ('<"••)•■/ A J 


or 


M„ 


/, r -m/eD\ + _L_ / y r -m/p0\ . — r -*/p0 + 

\ 2 / C0 2 \ C0 2 / pD (i) 


which reduces to 


— L- y + — 1 Y = — + C r m ^ p ^ 

V°2 M C0 2 S - ’ 


(A45) 


Evaluating equation ( A45 ) at r «* r^ and using equation (A40), 
yields 

J- Y 0 ♦J-V co •iUc.rf ’ 0 
”02 2b m co 2 C0 2b m 1 b 


Therefore the constant is 

2Y 0 Y co 

c . , °2b C0 2b A \ 1 

c i l Ter* u 


0 2 "C0 2 m j r m/pD 


157 


Substituting equation (A44) for the constant A in this equation, and 
substituting the resulting equation for the constant Cj and equa- 
tion (A44) into equation (A45), yields 




1 

XO. 


' CO , 


^(a*< 

m V 


'2$ 


* * C0 2s) 


♦ 



Using equation (A21), rb ■ 
( A46) can be written as 


'CO 


M 


2 J>. 


CO. 


W 2 * 



2S 

CO 



), equation 
2S 


(A46) 



(A47) 


Equation (A35) was solved exactly like (A34). Because the equa- 
tions and the boundary conditions are basically the same, only the 
result is presented: 



(A48) 


Equation (A36) was integrated, yielding 


2mc p T b 2XT b 
T ‘ H, 


XO 


XO 


dr 
dr " H 


pD 


m 

CO. 


CO. 


dY 


CO, 


M 


CO. 


CO. 

{ 

dr 


A" 


(A49) 


Writing this equation at the surface 


2mcT\ 2 aT k a 

E_k x k r 

H C0 S H C0 S dr 


P D 


M 


m 

CO, 


CO. 


'CO 


2$ 


M, 


CO. 


dY C0. 
r s dr' 


A" 




shows that 


A" 

or 

A" 


;mC 0 T b 

H co T$ 



2mcT. 

— E_b 

H co 5 



(V°B - "co^coj 4 ^C0 2s r s 

( 2H o 2 *o 2s - 2h co 2 *co 2S 4 *co 2s h co) 


The chemical formulae with the heat of reaction for the CO2 surface 
reaction and the gas phase reaction can be written as 
2C + 2 CO 2 «• 4C0, 2H C0 ^ 

2C0 + 0 2 2C0 2 , H cq 

Adding these equations yields 

2C ♦ 0 2 - 2C0, (h co - 2H C0 J 


Comparing this formula to the chemical formula with the heat of re- 
action for the O2 reaction at the surface, 


2C + 0 9 -► 2C0, H n , 
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shows that 


h o 2 “ ^ H C0 " 2H co 2 ^ 


Rewriting equation (A50) as 

A 1 


■ ■ ¥ x * * ft [v»« * (-*s + h 4 


(A51 ) 


and substituting equation (A51) into it, yields 


2mc_T h r s / 
A" m -jMu ♦A 

H C0 S CO \ 


2H.» *n + ^co 
°2 °2S C0 2S °2 




or 


A'' 


2mcJ‘ 


*^t, 


r S H 0, 


H co 5 H co 


( 2 %S + "Hs) 


(A52) 


Using the Schwab-Zeldovich approximation, D CQ . x/pc p * D, 
equation (A49) can be written as 

2mc. 


c p T b 2Dpc p T b dt m Y x p_D_ r 

£" T ‘ T r * ' "co 2 S "co/ 


dY 


CO, 


dr 


A" 


or 


(' m * p0r fe)*4r(* 


. 2 Vb 


dY 


mY C0« + p0r 


CO, 

dr 


A" 
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Dividing this equation through by pDr yields 


2c JV 


/ m . 1 / m y ®®2\ 

^ JBF T * #/* ^ JBr V C 0 j ♦ ~zr) 


H co V oDr 


Using the Integrating factor 
was Integrated, yielding 


. A" 
pDr 




(m/pDr)dr 


■ r" m ^ pD , this equation 


2c.J, 
“H 


^(> o0 )^(V n ' s0 )- / 


or 


2c_T, 
“H 


£ ("■“•*)• £ (V") 


A" 


’ D te) 


,-m/pO + 


which reduces to 


2c T. t C0 9 

b _ + £ 

H co m co 2 


- £ ♦ C.r"* D 
m * 


(A53) 


Evaluating equation (A53) at r = r b and using equation (A42), 
yields 


Iffilfe ♦ 3b . . 41 * c r" / < ,D 
h„ a : c 2 r b 


CO "CO, 


ir. 


Therefore the constant C 2 is 


C 2 ‘ 


^ + 3t 


H co M co, 


A" ] 1 

m / m/pO 


(A54) 


6 


Substituting equation (A5 2) for the constant A" into (A54), and sub 
stituting the resulting equation for C? along with equation (AS2) 
into (AS3), yields 
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Assuming the gas phase reaction 2C0 + 0^ > 2C0 2 is first order 
with respect to CO and 0^, the volume rate of consumption of CO can 
be written as 


^C0 “ k g C C0 C 0 2 eXP ^" E ^ RT ^ 


(A64) 


where k 

9 



R 

Letting 



« frequency or pre-exponential factor for the gas 
phase reaction CO + 1/2 0^ » C0 2 ; 

* molar concentration of carbon monoxide; 

* mo’ar concentration of oxygen; 

s activation energy for the gas phase reaction; 

CO ♦ 1/2 0 2 > C0 2 ; and 

* universal gas constant. 


(A65) 


and using the relationship 

pYf 

c i “ 

equation (A64) can be written as 


2 

*C° = k g Y C0 Y 0 2 exp( ‘ 6 l /x) 


(A66) 


Equation (A61), (A62), and (A63) can be substituted into equation 
(A66). The resulting equation can be substituted into the differen- 
tial equation for CO^, (A17), which then can be solved using a 
Runge-Kutta technique to obtain the C0 2 concentration as a function 
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of distance from the surface. However, because the CO^ concentra- 
tion should first increase with distance from the surface and then 
decrease, it was found that equation (A18) ft: ■ e CO concentration, 
which just decreases with distance from the surface, would be better 
to use in the Runge-Kutta technique. 

Equations (A61), ( A62 ) , and (A63) were solved so that _ 

f ' Y C0‘ r )* Y 0 = f ( Y C0’ r )’ and r= f ( Y C0’ r) ‘ result “ 
ing equations are 



(A69) 


The differential equation for the concentration of CO, equation (A18), 
can be written as 


m dY C0 eJD dY C0 
r ~5r~ " TT 


pO 


A 


dr 


jp = -M 


C0~C0 


or 





^co 

pD 
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Using 
_ m 

Y 

and rearranging, this equation becomes 



CO 

2 


M 


CO 




pD CO 


(1 - y) dY CO 
r dr 


(A70) 


Denoting equation (A67) and (A69) by Y Q = Y Q ( Y co* and 
T ~ r ^ Y C0’ r )» res P ect ’ vely, and substituting these into equation 
(A66), and substituting the result into equation (A70), yields 


d 2 V 


dr 


CO 

T 


k g OfyT Y C0 Y 0 ? (Y C0.r )exp [* 6 1 /T(y C0,r ) ] " (1 r ^ “dr^ < A71 ) 


If one of the variables Y fl Y CQ , Y co , or r is known, the 

others can be determined using equations (A61), (A62), and (A63). In 

this case, the surface temperature is considered to be known; there- 


fore, Y n , Y rn , and Y 
'2S LU 2S 


CO, 


can be determined. 


To determine these variables, first the consumption rates of 0^ 


and CO^ at the surface are written as 


*0 * k 0 ^0 

u 2 u 2 u 2 


and 

^C0 2 = k C0 2 C C0 2 


where 



= rate constant for the reaction, 2C + 0 2 
= rate constant for the reaction, C + C0 2 
= molar concentration of oxygen; and 


> 2C0; 
» 2C0; 




= molar concentration of carbon dioxide. 


'CO. 


Using the relationship 


pY 


1 


'1 M i ’ 


these equations become 


>Hf * _£_ k v 
°2S \ °2°2S 


(A72 


and 




C0 2S ' h C0 2 HHs 


(A73 


Substituting equations (A 72) and (A73) into (A21) yields 


V °2S Yc0 2S 

m * pM c r $ [2k 0 ^ ^ 


2 "CO, 


and, therefore, 

M,.r, 


Y = 


m 


rs 


pO 


2k 


2S 


’CO 


0, ft 


+ k 


2S 


2 " 0 . 


C0 2 M C0, 


(A7< 


? " 2 

Next equations (Abl) and (A63) are evaluated at the surface, yielding 





[ 26 (t s •'>♦«♦ »co a J • c r 

equations (A75) and (A76) can be written as 



Rewriting these equations as 



and subtracting, yields 
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Substituting equation (A79) into (A74), yields 


Y * 


M C r S 


(2V 


°2 C 1 


0 2 C 2 


M, 


CO, 


CO, 


'CO, 


r co 


2S 




(A80) 


Letting 


M C r S / 2k °2 C 1 

0 V% 2 C 2 ' 

and 

2k 0 2 M C r S ( C 1 

ICd 


M, 


CO, 


2M, 


CO, 

*C0, 


* C, 


■ c. 


equation (A80) can be written as 


C 3 Y C0 


+ C, 


2S 


(A81 ) 


Substituting equation (A81) into (A78) yields 


'CO 


* -e + C, 


2S 


ft) 


( C 3 Y C0 2S +C 4) 


(A82) 


which was solved by trial and error for Y rn . Once Y rn was 

LU 2S LU 2S 

obtained, equations (A61) and (A62) were evaluated at the surface and 


solved for Y 


2S 


and Y 


co s - 


The derivative of the mass fraction at the surface with respect 
to r can be obtained using equations (A24), (A26), and (A28). This 
will be illustrated using equation (A28). Substituting equations 
(A72) and (A73) into equation (A28) yields 
mY, 

+ 


_ _2_d_ ^co 


M co dr 


'CO 


r S M C0 


I , le_ k y 

h o 2 V«8 


X is k Y 

C0„ C0„«. 
C0 2 2 2S 




T^wiKiiMur mri 
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or 


dY 


CO 

dr 


m 


r $ pD CO $ D V M q 


M co (Ws , S Y ”a\ 

a ~^r) 


(A83) 


This equation was solved to obtain the derivative of the mass fraction 
of CO at the surface. 

With the mass fraction of CO, Y co , and its derivatave at the 
surface, { s * ^ e ^ n 9 known, equation (A71) can be written as 


F • k 


9 DMq w 


Wo ' V r)exp [- V ,(v co- r) ] • iL r 1£l “sr (A84) 


This equation was solved using a Runge-Kutta technique, by defining 
K1 * ArF 


K3 = Ar F 


4 Y ‘M 
l K * Y C0’ dr J 

(■ 

(n + 4 r Y Ar dY C0 . ArKl dY C0 . K2l 

l R + 2 ’ Y C0 + 2 ~3 F 8 * ~dr" + T 1 


K2 = ArF ( R + — Y + ^ ^ + Mi dY C0 . Kl] 
KZ ArF \ R + 2 . Y co + 2 + -5-, -^r* + -2) 


K4 - ArF + Ar, Y c() + Ar ^ ^ + K3^ 


Y C0 = Y C0 + Ar 


Ao 

dr 


+ (Kl + K2 + K3)/6 


3 


dY co dY co 

Tr “ "dr + ^ K1 + 2K2 * 2K3 + K4)/6 
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T he solution of equation (A84) gave values of the mass fraction of CO 
as a function of distance from the surface. The mass fractions of CO 
were then used in equations (Ab7), (A68), and (A69) to obtain the mass 
fractions of 0^ and CO^ and the dimensionless temperature, * , as a 
function of distance from the surface. 
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